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ABSTRACT 


Unmanned  Aircraft  Systems  (UASs)  are  critical  for  future  combat  effectiveness. 
Military  planners  from  all  branches  of  the  Department  of  Defense  now  recognize  the 
value  that  real  time  intelligence  and  surveillance  from  UASs  provides  the  battlefield 
commander.  The  Operations  Analysis  Division  of  the  Marine  Corps  Combat 
Development  Command  is  currently  conducting  an  Overarching  Unmanned  Aircraft 
Systems  study  to  determine  future  force  requirements.  Current  analysis  is  conducted 
through  the  use  of  the  Assignment  Scheduling  Capability  for  Unmanned  Air  Vehicles 
(ASC-U)  and  several  specially  designed  heuristics.  The  Unmanned  Aircraft  System 
Scheduling  Tool  (UAS-ST)  combines  these  capabilities  into  one  model  and  addresses 
several  issues  associated  with  ASC-U.  UAS-ST  allows  the  user  to  control  all  aspects  of 
the  UAS,  define  a  scenario,  and  then  generates  a  flight  schedule  over  a  known  time 
horizon  based  on  those  inputs.  All  missions  are  assigned  a  user  defined  value  and  the 
total  schedule  value  is  reported.  The  user  can  then  quickly  change  a  parameter  of  the 
UAS,  re-solve  the  model,  and  see  the  impact  their  proposed  change  has  on  the  overall 
value  of  the  schedule  attained.  Therefore,  UAS-ST  is  a  tool  for  analyzing  the  value  of 
future  changes  in  UAS  structure. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


THESIS  DISCLAIMER 


The  reader  is  cautioned  that  computer  programs  developed  in  this  research  may 
not  have  been  exercised  for  all  cases  of  interest.  While  every  effort  has  been  made, 
within  the  time  available,  to  ensure  that  the  programs  are  free  of  computational  and  logic 
errors,  they  cannot  be  considered  validated.  Any  application  of  these  programs  without 
additional  verification  is  at  the  risk  of  the  user. 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


TABLE  OF  CONTENTS 


I.  INTRODUCTION . 1 

A.  PURPOSE  AND  OVERVIEW . 1 

B.  BACKGROUND . 2 

1.  Problem  Statement . 2 

2.  The  Marine  Corps  Three  Tier  UAS  Family  of  Systems . 3 

3.  Current  UAS  Structure . 5 

4.  Current  UAS  Missions  Types . 6 

C.  SCOPE  AND  LIMITATIONS . 7 

D.  THESIS  ORGANIZATION . 7 

II.  ASSIGNMENT  SCHEDULING  CAPABILITY  FOR  UAS . 9 

A.  INTRODUCTION . 9 

B.  BASICS  OF  ASC-U . 9 

1.  Model  Description . 9 

2.  ASC-U  Implementation . 10 

3.  Model  Capabilities  and  Limitations . 10 

C.  LITERATURE  REVIEW . 12 

III.  OPTIMIZATION  MODEL  FOR  THE  SCHEDULING  OF  UNMANNED 

AIRCRAFT  SYSTEMS . 15 

A.  INTRODUCTION . 15 

B.  AN  INTEGER  PROGRAM  TO  OPTIMIZE  UAS  SCHEDULING . 15 

1.  Indices . 15 

2.  Given  Data  [units] . 16 

3.  Decision  Variables . 17 

4.  Formulation . 18 

5.  Discussion . 19 

6.  Route  Enumeration  and  Data  Development . 19 

IV.  COMPUTATIONAL  RESULTS . 23 

A.  INITIAL  TE  ST  SCENARIO . 23 

B.  UAS-ST  INITIAL  RESULTS  AND  ANALYSIS . 25 

C.  COMPARING  UAS-ST  TO  ASC-U . 29 

V.  CONCLUSIONS  AND  NEW  OPPORTUNITIES . 33 

A.  SUMMARY . 33 

B.  OPERATIONAL  INTRODUCTION . 33 

C.  FUTURE  DEVELOPMENT . 34 

LIST  OF  REFERENCES . 35 

APPENDIX  A  TEST  SCENARIO . 37 

APPENDIX  B  VBA  CODE . 41 

APPENDIX  C  GAMS  CODE . 65 

INITIAL  DISTRIBUTION  LIST . 75 


IX 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


x 


LIST  OF  FIGURES 


Figure  1.  UAS  Tier  Structure . 3 

Figure  2.  Projected  Future  FoS  Tiers  Structure . 5 

Figure  3.  Plot  of  OUAS  Scenario  Taken  from  OUAS . 24 

Figure  4.  Example  of  Payload  Worksheet . 27 

Figure  5.  OUAS  Study  Future  Tier  Perfonnance  Parameters . 27 

Figure  6.  Sample  of  User  Interface  (Dashboard) . 28 

Figure  7.  Partial  Display  of  Optimal  UAS-ST  Schedule . 29 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


LIST  OF  TABLES 


Table  1 .  General  Description  of  Mission  Types . 6 

Table  2.  Number  of  UAS  Assets  Available  in  Generic  Scenario . 23 

Table  3.  Input  Worksheet  for  Mission  Requests . 25 

Table  4.  Assigned  Mission  Values  by  Tier . 26 


xiii 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


xiv 


EXECUTIVE  SUMMARY 


Unmanned  Aircraft  Systems  (UASs)  procurement  is  vital  to  the  United  States 
Marine  Corps'  (USMC)  combat  effectiveness  in  the  near  future.  UASs  are  used  for 
collecting  intelligence,  surveillance,  and  targeting  infonnation.  They  accomplish  these 
missions  at  a  much  lower  overall  risk  than  conventional,  manned  aircraft.  Military 
planners  from  all  branches  of  the  Department  of  Defense  (DoD)  recognize  the  value  that 
real  time  intelligence  and  surveillance  from  UASs  provides  the  battlefield  commander. 
In  an  ongoing  operation,  mission  requests  from  units  in  theater  typically  far  exceed  the 
capacity  of  available  UAS  assets.  Demand  for  UAS  missions  is  increasing  as  the 
capability  of  these  platforms  expands.  Individual  branches  of  DoD  are  scrambling  to 
acquire  the  UASs  needed  to  support  the  requirements  of  currently  deployed  units  in 
combat  operations.  As  a  result,  interoperability  and  compatibility  is  a  major  concern  with 
today’s  current  family  of  UASs. 

Due  to  operational  requirements  in  Iraq  and  Afghanistan,  DoD  is  focusing 
procurement  strategy  away  from  force  transformation  and  future  generation  weapons  to 
more  immediate  concerns.  Acquisition  of  UASs  needed  to  meet  operational  requirements 
is  the  agency’s  highest  priority.  As  a  result,  large  procurement  budgets  exist  to  fill  the 
supply  shortage  and  meet  future  operational  requirements.  Currently  DoD  is  evaluating  a 
large  number  of  alternative  UAS  programs.  This  selection  process  requires  an  objective 
analysis  of  each  alternative.  The  model  in  this  thesis  allows  the  USMC  to  accomplish 
their  analysis  by  clearly  depicting  the  impact  of  changes  in  system  capabilities  on  a  daily 
flight  schedule.  Application  of  the  model  is  not  limited  to  UASs.  The  model  is  capable 
of  analyzing  any  asset  which  must  be  scheduled  in  response  to  user  demands.  Therefore, 
this  model  has  future  application  in  a  multitude  of  other  programs. 

The  current  model  in  use  for  the  Marine  Corps  Combat  Development  Command 
(MCCDC)  Overarching  Unmanned  Aircraft  Systems  (OUAS)  study  is  the  Assignment 
Scheduling  Capability  for  UAVs  (ASC-U)  model.  The  design  agency  for  ASC-U  is  the 
U.S.  Army  Training  and  Doctrine  Command  Analysis  Center  (TRAC).  ASC-U  is 
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designed  to  support  the  development  of  an  effective  UAS  force  structure.  It  is  a 
spreadsheet-based  decision  support  tool  primarily  for  allocation  and  scheduling  of  assets. 
ASC-U  addresses  complexities  in  military  operations  and  scheduling  of  multiple  moving 
platforms.  ASC-U  accepts  user  parameters  that  define  a  scenario  then  seeks  to  provide  a 
feasible  schedule  for  available  UASs.  ASC-U  is  the  first  model  that  enables  analysts  to 
address  this  scheduling  problem  effectively.  ASC-U  combines  both  optimization  and 
simulation  to  produce  a  tool  with  unique  capabilities. 

This  thesis  develops  an  integer  linear  programming  model,  UAS-ST,  for 
scheduling  UASs.  UAS-ST  allows  the  user  to  define  all  elements,  both  operational  and 
performance,  of  the  UAS  via  an  Excel  spreadsheet.  A  schedule  generator  written  in 
visual  basic  for  this  application  then  takes  these  elements  and  generates  a  user-defined 
number  of  individual  schedules.  This  schedule  generation  is  done  for  every  UAV  (UAV) 
that  is  included  in  the  scenario.  Once  the  predefined  number  of  schedules  is  generated, 
UAS-ST  creates  data  files  for  a  General  Algebraic  Modeling  System  (GAMS)  model. 
GAMS,  through  the  use  of  CPLEX,  finds  a  near  optimal  combination  of  individual 
schedules  to  produce  a  complete  schedule  for  a  user  designated  time  period. 

The  initial  test  scenario  replicates  a  generic  scenario  given  in  the  OUAS  study. 
This  replication  provides  the  simplest  and  most  direct  comparison  of  results  with  ASC-U. 
The  scenario  time  period  is  twenty  four  hours  divided  into  ninety  six  intervals  of  fifteen 
minutes  each.  This  scenario  does  not  refer  to  a  specific  country,  but  is  designed  solely  to 
provide  the  framework  for  determining  the  value  of  changes  in  UAS  structure  on  a  given 
set  of  mission  requests.  UAS  force  structure  for  the  initial  scenario  represents  a  small  but 
realistic  composition  of  systems.  Tier  I  is  comprised  of  three  separate  units  with  four 
UAVs  and  one  Ground  Control  Station  (GCS)  per  unit.  Tier  II  is  comprised  of  three 
separate  units  with  one  UAV  and  one  GCS  per  unit.  Tier  III  is  comprised  of  a  single  unit 
with  two  UAVs  and  two  GCSs. 

UAS-ST  is  significantly  different  from  ASC-U.  ASC-U  utilizes  an  Excel 
spreadsheet  and  pulls  the  data  into  an  Access  database.  UAS-ST  applies  an  Excel 
spreadsheet,  which  then  uses  Visual  Basic  to  transfer  the  data  to  GAMS.  GAMS  then 
applies  the  CPLEX  solver  to  quickly  optimize  the  schedule  for  all  mission  requests. 


Features  that  require  supplemental  heuristics  in  ASC-U  are  incorporated  directly  into 
UAS-ST,  eliminating  the  need  for  further  processing.  This  produces  a  much  quicker  and 
efficient  analysis  of  alternatives.  CPLEX  is  highly  efficient  in  preprocessing  feasible 
solutions  and  reduces  the  run  time  to  a  matter  of  minutes.  Overall,  UAS-ST  provides  an 
efficient  update  to  the  model  presented  in  ASC-U. 

UASs  are  critical  to  our  nation’s  and  the  USMC's  future  military  combat 
effectiveness.  All  branches  of  DoD  recognize  the  need  to  develop  an  integrated  network 
of  UASs.  To  answer  this  need,  the  Marine  Corps  is  developing  the  concept  of  the  UAS 
Family  of  Systems  (FoS).  The  FoS  calls  for  a  three  tier  structure  of  UASs  with 
overlapping  capabilities.  Currently,  the  MCCDC  Operation  Analysis  Division  is 
conducting  an  OUAS  study.  This  thesis  is  a  direct  contribution  to  that  study,  and  UAS- 
ST,  provides  a  planning  tool  for  development  of  future  UAS  structure.  UAS-ST  allows 
the  user  to  control  all  aspects  of  the  UAS,  define  a  scenario,  and  then  generates  a  flight 
schedule  over  a  known  time  horizon  based  on  those  inputs.  All  missions  are  assigned  a 
user-defined  value  and  the  total  schedule  value  is  reported.  The  user  can  then  quickly 
change  an  operational  or  perfonnance  parameter  of  the  UAS,  re-solve  the  model,  and  see 
the  impact  on  the  overall  value  of  the  schedule.  UAS-ST  is  a  tool  for  analyzing  the  value 
of  future  changes  in  UAS  structure. 
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I.  INTRODUCTION 


A.  PURPOSE  AND  OVERVIEW 

Unmanned  Aircraft  System  (UAS)  procurement  is  vital  to  the  United  States 
Marine  Corp’s  (USMC)  combat  effectiveness  in  the  near  future.  UASs  are  used  for 
collection  of  intelligence,  surveillance,  and  targeting  infonnation.  They  accomplish  these 
missions  at  a  much  lower  risk  than  conventional,  manned  aircraft.  Military  planners  from 
all  branches  of  the  Department  of  Defense  (DoD)  recognize  the  value  that  real  time 
intelligence  and  surveillance  from  UASs  provides  the  battlefield  commander.  In  an 
ongoing  operation,  mission  requests  from  units  in  theater  typically  far  exceed  the 
capacity  of  available  UAS  assets.  Future  demand  for  UAS  missions  will  increase  as  the 
capability  of  these  versatile  platfonns  expands.  Individual  branches  of  DoD  are 
scrambling  to  acquire  the  UASs  needed  to  support  the  requirements  of  currently  deployed 
units  in  combat  operations.  [IHT,  2008]  As  a  result,  interoperability  and  compatibility  is 
a  major  concern  with  today’s  current  family  of  UASs. 

To  address  these  issues  in  future  systems,  the  USMC  is  developing  the  UAS 
Family  of  Systems  (FoS).  Each  system  consists  of  an  Unmanned  Air  Vehicle  (UAV), 
Ground  Control  Station  (GCS),  Launch  and  Recovery  Station  (LRS)  and  a  combination 
of  various  sensors  and  payload  components.  The  components  of  the  FoS  have  several 
complementary  capabilities  which  overlap  in  certain  mission  areas.  The  intent  is  to 
create  a  mix  of  several  UASs  able  to  support  various  units  of  different  sizes  and  levels  of 
operation. 

To  assess  the  effectiveness  of  a  combination  of  UASs,  quantitative  models  must 

be  applied  to  provide  a  reasonably  accurate  measurement  of  capability  and  some 

guidance  for  efficient  employment.  To  this  end,  this  thesis  develops  a  UAS  planning  and 

decision  support  tool  that  takes  as  input  a  planning  horizon,  a  fleet  of  UAVs  and  their 

individual  operating  limits,  a  list  of  available  payloads  and  a  list  of  mission  requests.  The 

model  then  provides  as  output  an  operational  schedule  for  each  individual  vehicle 

indicating  the  missions  to  be  accomplished  in  a  specified  time  horizon  and  the  payloads 
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required.  This  schedule  provides  analysts  the  ability  to  quickly  determine  the  impact  of 
changes  in  any  system  parameter  on  the  overall  value  of  missions  accomplished. 

B.  BACKGROUND 

1.  Problem  Statement 

Due  to  operational  requirements  in  Iraq  and  Afghanistan,  DoD  is  focusing 
procurement  strategy  away  from  force  transformation  and  future  generation  weapons  to 
more  immediate  concerns.  Acquisition  of  the  UASs  needed  to  meet  operational 
requirements  is  the  agency’s  highest  priority.  [IHT,  2008]  As  a  result,  large  procurement 
budgets  exist  to  fill  the  supply  shortage  and  meet  future  operational  requirements. 
Currently  DoD  is  evaluating  a  large  number  of  alternative  UAS  programs.  The  selection 
process  requires  an  objective  analysis  of  each  alternative.  This  thesis  develops  a  model 
to  help  the  USMC  with  this  analysis.  The  model  in  this  thesis  allows  the  USMC  to 
accomplish  their  analysis  by  clearly  depicting  the  impact  of  changes  in  system 
capabilities  on  a  daily  flight  schedule.  Application  of  the  model  is  not  limited  to  UASs. 
The  model  is  capable  of  analyzing  any  asset  which  must  be  scheduled  in  response  to  user 
demands.  Therefore,  this  model  has  future  application  in  a  multitude  of  other  programs. 

In  recent  comments,  Defense  Secretary  Robert  Gates,  addresses  wasteful  or 
inefficient  UAS  procurement  programs  in  a  speech  given  at  the  Air  Force’s  Air 
University  at  Maxwell  Air  Force  Base: 

Because  people  were  stuck  in  old  ways  of  doing  business,  it's  been  like 
pulling  teeth.  While  we've  doubled  this  capability  in  recent  months,  it  is 
still  not  good  enough.  [IHT,  2008] 

The  fact  that  the  Defense  Secretary  is  unusually  blunt  in  his  criticism  of  current  program 
development,  should  serve  as  proof  of  the  pressure  to  field  future  UASs. 
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2. 


The  Marine  Corps  Three  Tier  UAS  Family  of  Systems 


USMC  Three-Tier  UAS  Strategy 

Ttor  III:  VUAB 
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•  Shadow  aaOCOuntt) 
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Figure  1.  UAS  Tier  Structure 


The  FoS  divides  the  various  UASs  into  three  separate  tiers.  The  tier  assignments 
provide  each  level  of  the  Marine  Air  Ground  Task  Force  (MAGTF)  an  organic, 
interoperable,  integrated  and  tailored  capability  that  raises  the  situational  awareness  of 
the  combat  unit  commander  through  a  common  communication  network.  The  current 
operational  conditions  in  Operation  Enduring  Freedom  and  Operation  Iraqi  Freedom 
demonstrate  the  importance  of  UAS  operations  and  the  need  for  commanders  at  all  levels 
to  maintain  control  of  their  respective  battle  space.  Therefore,  tiers  are  defined  by  the 
area  of  interest  and  operational  level  of  the  supporting  unit,  such  as  a  company,  battalion, 
or  regiment.  There  is  some  operational  overlap  in  tier  capabilities;  this  is  intentional  and 
should  be  considered  beneficial  to  mission  accomplishment. 
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The  three  tiers  are  described  in  [MCCDC,  2005]  as  follows: 

■  Tier  I  -  Short  Range  UAS 

•  Operationally  Supports  Battalion,  Company,  and  Platoon 

•  Performs  Reconnaissance  and  Surveillance  Missions 

•  Current  System:  Dragon  Eye,  Future:  Raven-B 

•  Endurance  of  1.5  hours,  Combat  Radius  of  5  miles 

•  Speed:  Less  than  20  knots 

•  Capacity  for  Single  Payload 


■  Tier  II  -  Division  or  Regimental  UAS 

•  Operationally  Supports  Division,  Marine  Expeditionary  Unit 
(MEU),  Regiment  and  Battalion 

•  Air  vehicle:  Persistent,  low  cost,  durable,  low  observable,  easily 
transported,  shipboard  compatible,  target  acquisition  capable, 
heavy  fuel  engine 

•  Current:  Scan  Eagle  Contract 

•  Endurance  of  15  hrs,  Combat  Radius  of  50  miles 

•  Speed:  60  knots 

•  Capacity  for  Two  Payloads 


-  Tier  III  -  UAS 

•  Operationally  supports  Marine  Expeditionary  Force,  Marine 
Expeditionary  Brigade,  MEU,  Division  and  Regiment 

•  Reconnaissance  Surveillance  and  Target  Acquisition,  Electro- 
Optical  /Infrared  Imagery  (EO/IR) 

•  Current  System:  Pioneer,  Future  Concept  -  Vertical  Takeoff  UAV 

•  Air  vehicle  (Proposed;  expeditionary,  sea  based,  Vertical  Takeoff, 
large  sensor  payload,  multi-mission  capable,  weaponization) 

•  Sensor  payload  (Proposed:  EO/IR/laser  designation, 

communication  relay,  Signals  Intelligence,  Electronic  Warfare) 

•  Endurance  of  8  hours,  Combat  Radius  of  300  miles 

•  Speed:  200  knots 

•  Capacity  for  three  payloads 
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FoS  Tiers 


Today 

Tier  I 

•  Dragon  Eye 


2020 

Tier  I 

•  Joint  Manpackable  UAS  (Raven  B) 

•  4  systems  per  inf  battalion 

•  System  =  1  GCSw/  3  UAVs 


Tier  II 

•  ISR  services  contract  (Scan  Eagle) 


Added 

ca£abilities 


> 


Tier  II 

•  TBD  endurance  system 

•  9  systems  per  division 

•  System  =  1  GCS  w/3  UAVs 


Tier  III 

•  Pioneer 

•  2  VMUs 

•  1  system  perVMU 

•  System  =  1  GCSw/ 5  UAVs 


Tier  III 

•  VUAS 

•  4 VMUs 

•  3  systems  per  VMU 

•  System  =  2  GCSw/  4  UAVs 


Figure  2.  Projected  Future  FoS  Tiers  Structure 


3.  Current  UAS  Structure 

There  are  several  specialized  terms  associated  with  a  UAS  that  this  thesis  uses 
repeatedly.  To  avoid  any  confusion,  the  following  definitions  apply: 

•  UAV  (Unmanned  Aerial  Vehicle):  A  UAV  is  an  unpiloted  aircraft.  UAVs 
can  be  remote  controlled  or  fly  an  automated  route  based  on  pre-programmed 
information. 

•  UAS  (Unmanned  Aircraft  System):  UAS  is  the  current  term  introduced  by 
DoD  and  accepted  by  the  Federal  Aviation  Administration  to  replace  the  tenn 
UAV.  A  UAS  consists  of  not  only  the  unmanned  aircraft,  but  also  the  data 
link  system,  the  launch  and  recovery  station,  and  all  maintenance  and  support 
equipment. 

•  GCS  (Ground  Control  Station):  A  GCS  is  a  land  or  sea-based  system  that 
allows  an  operator  to  control  an  unmanned  aircraft.  A  GCS  may  control  both 
"active"  and  "passive"  missions.  An  active  mission  requires  the  operator  to 
monitor  the  UAV,  while  a  passive  mission  requires  no  interaction  from  the 
operator. 
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•  LRS  (Launch  and  Recovery  Site):  Used  to  control  an  unmanned  aircraft 
during  the  initial  and  terminal  phases  of  flight.  The  LRS  can  also  function 
as  a  GCS  when  not  recovering  UAVs. 

•  Payload:  Future  system  requirements  call  for  a  system  of  “payloads,” 
each  of  which  provides  a  custom  capability.  These  payloads  will  be 
uniform  in  size  and  connectivity  to  allow  for  rapid  configuration  of 
mission  specific  profiles. 

•  Mission  Package:  Consists  of  the  combination  of  payloads  loaded  on  an 
individual  UAV  which  determines  the  mission  capability  of  the  UAV. 
Required  payload  capacity  is  a  key  element  of  the  current  study. 

4.  Current  UAS  Missions  Types 

Following  is  a  table  of  terms  associated  with  current  UAS  missions.  These 
missions  represent  general  categories  and  are  derived  from  the  Overarching  Unmanned 
Aircraft  System  (OUAS)  study  [OUAS,  2007]. 


Missions 

Abbreviation 

Description 

Reconnaissance,  Surveillance,  Target 
Acquisition  (RSTA)  (EO/IR) 

RSTA-EO/IR 

ISR  asset  for  routine  day  and  night  time 
operations 

Reconnaissance,  Surveillance,  Target 
Acquisition  (Synthetic  aperture  radar) 

RSTA-SAR 

ISR  asset  for  dense  vegetation  and  poor 
weather  conditions 

Signals  Intelligence 

SIGINT 

Sensor(s)  designed  for  passive  collection  of 
signals 

Air  Vehicle  Communications  Link  Relay 

AVCLR 

Relays  information  and  instructions  through 
one  UAV  to  another 

Communications  Relay 

CR 

Relays  voice  and  data  between  ground  points 

Strike 

STK 

Weapons  enabled  kinetic  destruction  of  a  time 
sensitive  target 

Electronic  Warfare 

EW 

Active  denial  of  radio  frequency 

Table  1.  General  Description  of  Mission  Types 
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C.  SCOPE  AND  LIMITATIONS 

The  intent  of  this  thesis  is  to  provide  a  model  for  use  as  an  analytical  tool  in 
determining  the  future  force  structure  for  UASs.  The  model  allows  the  user  to  quickly 
change  both  operational  and  perfonnance  parameters  for  each  UAS  tier.  Once 
parameters  are  set,  the  model  quickly  generates  a  near-optimal  schedule.  This  rapid 
schedule  generation  allows  an  analyst  to  see  the  impact  of  changing  UAS  capabilities  on 
a  given  daily  schedule. 

D.  THESIS  ORGANIZATION 

Chapter  II  provides  a  discussion  of  the  Assignment  Scheduling  Capability  for 
UAVs  (ASC-U)  the  current  UAS  evaluation  tool  in  use  by  the  Operation  Analysis 
Division  (OAD),  MCCDC.  Chapter  III  describes  the  optimization  model  and  the  stack 
based  enumeration  heuristic  used  to  solve  it.  Chapter  IV  provides  a  detailed  analysis  of 
the  ASC-U  results  and  compares  it  to  the  current  model.  Chapter  V  is  devoted  to 
conclusions  and  recommendations  for  future  research. 
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II.  ASSIGNMENT  SCHEDULING  CAPABILITY  FOR  UAS 


A.  INTRODUCTION 

The  current  model  in  use  for  the  OUAS  study  is  the  ASC-U  model.  The  agency 
responsible  for  the  design  of  ASC-U  is  the  U.S.  Army  Training  and  Doctrine  Command 
Analysis  Center  (TRAC).  ASC-U  is  designed  to  support  the  development  of  an  effective 
UAS  force  structure.  This  chapter  describes  the  basic  structure  of  the  ASC-U  model.  It 
begins  with  a  description  of  the  inputs  accepted  by  the  model  then  describes  its 
implementation  in  the  study.  The  chapter  then  discusses  limitations  of  ASC-U  which  are 
addressed  by  this  thesis.  The  majority  of  the  information  given  here  is  summarized  from 
the  OUAS  report  [OUAS,  2007].  This  chapter  also  includes  a  review  of  current 
literature. 

B.  BASICS  OF  ASC-U 

1.  Model  Description 

ASC-U  supports  the  development  of  an  effective  UAS  force  structure.  It  is  a 
spreadsheet-based  decision  support  tool  used  primarily  for  allocation  and  scheduling. 
ASC-U  addresses  the  complexities  in  military  operations  and  scheduling  of  multiple 
moving  platforms.  ASC-U  accepts  user  parameters  that  define  a  scenario  and  then  seeks 
to  provide  a  feasible  schedule  for  available  UASs.  As  stated  in  the  ASC-U  Analyst 
Manual: 

ASC-U  provides  a  solution  to  the  following  problem:  Given  a  scenario 
that  specifies  the  number  of  each  type  of  UAV,  initial  UAV  locations,  and 
UAV  performance  characteristics,  detennine  the  number  of  missions  that 
can  successfully  be  completed  and  the  schedule  for  each  UAV.  The 
solution  must  consider  GCS  locations  and  capacities,  remote  viewing 
tenninal  requirements,  and  communication  platfonn  footprint  and 
capacities.  [Ahner,  2006] 
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ASC-U  is  the  first  model  that  enables  analysts  to  address  this  scheduling  problem 
effectively.  It  combines  both  optimization  and  simulation  to  produce  a  tool  with  unique 
capabilities. 

2.  ASC-U  Implementation 

The  ASC-U  model  works  by  allowing  the  user  to  design  a  scenario  which  consists 
of  a  set  of  mission  requests.  A  mission  request  consists  of  a  specific,  required  UAS 
capability  at  a  specific  geographic  location,  for  a  set  amount  of  time.  The  location  of  the 
mission  remains  fixed  once  it  is  assigned.  The  allocation  tool  uses  the  UAS  capability 
data  and  mission  data  it  is  given  to  create  a  feasible  schedule  which  will  accomplish  as 
many  mission  requests  as  possible.  ASC-U  is  deterministic.  For  a  given  input,  it  will 
always  produce  the  same  schedule  with  the  same  measures  of  performance.  Specifically, 
ASC-U  uses  a  detenninistic  algorithm  to  optimize  over  a  given  finite  time  horizon  to 
obtain  near-optimal  UAS  mission  area  assignments. 

ASC-U  allows  the  user  to  define  several  different  input  parameters.  Mission 
requirements  are  the  most  important  parameter.  Mission  requirements  consist  of  a 
coordinate  location,  payload  requirement,  length  of  mission  and  mission  priority. 
Mission  priority  functions  as  a  selection  criterion.  A  user  may  establish  mission 
precedence  by  assigning  a  higher  value  to  a  specific  mission  type.  UAS  parameters  that 
may  be  entered  by  the  user  are  payload  details,  GCS  attributes,  and  UAV  attributes.  GCS 
attributes  consist  of  coordinate  location,  control  limits,  and  unit  assignment.  UAV 
attributes  include  speed,  operating  time,  combat  radius,  launch  site,  and  total  time 
available. 


3.  Model  Capabilities  and  Limitations 

ASC-U  is  the  first  attempt  to  create  a  specific  tool  that  involves  all  aspects  of  the 
UAS  family  of  systems.  As  part  of  its  support  program,  TRAC-Monterey  publishes  an 
ASC-U  user’s  manual.  For  its  use  in  the  OUAS  study,  some  settings  are  different  than 
the  manual’s  recommended  settings.  Different  settings  are  required  due  to  the  small  size 

of  the  original  test  scenario  chosen.  ASC-U  is  designed  to  model  several  thousands  of 
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missions  over  a  long  time  period.  The  scenarios  used  in  the  OUAS  study  are  much 
shorter  in  duration  and  involve  only  a  few  hundred  missions. 

The  mission  requirements  are  the  most  important  input  data  because  they  have  the 
most  direct  effect  on  the  schedule  generated.  ASC-U  attempts  to  complete  as  many 
mission  hours  as  possible.  The  objective  function  has  a  significant  drawback.  It  leads  to 
preferential  assignment  of  missions  that  are  close  to  the  launch  and  recovery  site  because 
they  allow  for  more  follow-on  missions  to  be  accomplished.  The  model  schedules  as 
many  close-in  missions  as  possible  and  foregoes  farther  outlying  missions.  Therefore, 
mission  priority,  a  specific  input,  allows  the  user  to  set  a  precedence  level.  However,  in 
some  instances  the  user  is  forced  to  artificially  inflate  the  value  of  a  mission  to  ensure  it 
is  scheduled. 

All  inputs  are  usually  entered  in  Excel  spreadsheets,  and  then  read  into  Access. 
ASC-U  can  take  input  from  either  Excel  or  Access.  The  Excel  inputs  are  placed  in  an 
Access  data  base  as  the  model  begins  its  run.  The  output  is  also  stored  in  an  Access  data 
base.  Several  tables  are  produced  in  the  output.  The  most  important  table  for  the  OUAS 
study  is  mission  coverage.  Mission  coverage  is  broken  down  into  requesting  unit, 
mission  type  and  UAV  type.  In  ASC-U  every  possible  payload  combination  is 
enumerated  as  a  mission  package.  Mission  packages  are  then  assigned  to  UAVs. 
Mission  package  usage  is  recorded;  therefore  individual  payload  usage  is  not  available. 

The  objective  of  the  first  OUAS  study  is  to  provide  as  many  hours  of  UAS 
support  as  possible.  For  ASC-U,  the  two  key  factors  are  the  optimization  interval  and  the 
time  horizon  for  scheduling  the  UAV.  The  optimization  interval  controls  how  often 
ASC-U  runs  its  routine  to  assign  missions.  The  time  horizon  controls  how  far  forward  in 
time  ASC-U  will  look  to  assign  the  UAV  to  a  mission.  The  optimization  interval  is 
crucial  because  if  the  interval  is  set  too  long  then  the  UAS  will  remain  idle  instead  of 
performing  another  mission.  If  the  interval  is  set  too  low  the  model  has  difficulties  and 
exhibits  odd  behavior.  An  interval  of  12  minutes  is  best  for  the  scenarios  used  in  the 
OUAS  study  [OUAS,  2007].  Time  horizon  is  critical  because  ASC-U  will  only  schedule 
a  UAV  once  during  the  given  time  period.  Therefore,  if  the  horizon  is  set  too  long,  ASC- 

U  can  see  a  high  value  mission  at  the  end  of  the  period  and  keep  a  UAV  idle  the  entire 
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period  waiting  for  that  mission.  A  UAV  can  complete  another  mission  and  still  complete 
its  priority  mission.  On  the  other  hand,  if  the  time  horizon  is  set  too  low  then  ASC-U 
misses  high  priority  missions  because  it  cannot  see  far  enough  out.  The  final  horizon  for 
the  OUAS  study  is  6  hours  for  Tier  II  and  4  hours  for  Tier  III. 

Because  of  these  shortcomings,  several  workarounds  and  new  heuristics  are  used 
to  enable  ASC-U  to  find  schedules  that  exhibit  required  behaviors.  The  SUPER 
MISSION  ability  allows  a  UAV  to  fulfill  multiple  missions  if  they  are  within  its  operating 
range.  Three  optional  heuristics  are  also  used  in  the  OUAS  study.  They  are  implemented 
in  this  order: 

EARLY  RETURN  (Go  Home  vs.  Stay)  -  When  a  UAV  completes  an  assigned 
mission,  a  value  of  TRUE  allows  the  UAV  to  return  to  base  if  it  can  do  so  and  still  make 
it  back  on  station  in  time  for  its  next  mission.  A  value  of  FALSE  forces  the  UAV  to  fly 
until  its  operating  limit  is  met.  For  the  OUAS  study,  this  is  set  to  FALSE. 

SECONDARY  AREAS  (Go  Get  More  Value  From  Another  Mission)  -  If  a  UAV  is 
scheduled  for  more  than  one  mission  in  the  same  location  with  a  time  gap  in  between 
them,  a  value  of  TRUE  allows  the  UAV  to  perform  another  mission  in  between  as  long  as 
it  is  available  at  the  start  of  its  previously  scheduled  mission.  A  value  of  FALSE  will 
force  it  to  remain  on  station  until  the  start  of  its  next  mission.  For  the  OUAS  study,  this 
was  set  to  TRUE. 

APPENDED  AREAS  (Done,  Is  There  Another  Mission?)  -  If  a  UAV  has 
completed  its  assigned  missions,  a  value  of  TRUE  allows  it  to  use  its  remaining  time  to 
find  another  mission.  A  value  of  FALSE  will  not  allow  additional  missions.  For  the 
OUAS  study,  this  was  set  to  TRUE. 

C.  LITERATURE  REVIEW 

Many  recent  studies  attempt  to  shape  some  aspect  of  the  future  design  of  UAS 
force  structure.  An  unintended  consequence  is  that  most  of  these  studies  focus  only  on 
some  specific  technical  portion  of  the  overall  problem.  The  Deputy  Commandant, 
Aviation  and  the  Deputy  Commandant  Combat  Development  and  Integration  are 
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sponsoring  the  OUAS  study.  Their  intent  is  to  analyze  the  future  USMC  UAS  force 
structure  to  detennine  how  to  best  meet  the  needs  of  the  MAGTF.  The  initial  phase  of 
OUAS  is  being  conducted  by  OAD  MCCDC.  The  model  given  in  this  thesis  is  intended 
for  use  in  the  follow  on  analysis  of  OUAS.  Therefore,  the  initial  report  by  OAD  is 
critical  to  this  thesis,  as  it  largely  detennined  the  requirements  for  the  model. 

Tutton  [2003]  deals  with  the  optimal  placement  of  a  unit’s  sensing  assets.  He 
presents  a  methodology  for  finding  the  most  beneficial  mix  and  allocation  strategy  for  an 
individual  unit’s  sensors  for  a  given  threat  scenario.  Doll  [2004]  takes  the  model 
developed  by  Tutton  [2003]  and  translates  it  into  a  programming  language  for  easier 
simulation.  She  refines  many  of  the  constraints  in  the  original  model  to  make  for  a  much 
more  realistic  simulation.  Finally,  Zacherl  [2006]  deals  specifically  with  reactive  aircraft 
scheduling.  The  thesis  develops  a  model  which  reviews  a  current  air  tasking  order  and 
then  rapidly  reassigns  aircraft  to  new  targets  as  they  become  available. 

Finally  a  large  number  of  commercial  information  sources  address  current  topics 
in  UAS  development  and  model  optimization.  When  possible,  these  sources  are  used  for 
the  sake  of  currency  or  accuracy. 
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III.  OPTIMIZATION  MODEL  FOR  THE  SCHEDULING  OF 
UNMANNED  AIRCRAFT  SYSTEMS 


A.  INTRODUCTION 

This  chapter  develops  the  mathematical  programming  model,  Unmanned  Aircraft 
System  Scheduling  Tool  (UAS-ST).  UAS-ST  applies  an  optimization  based  approach  to 
analyze  the  best  mix  of  future  UAS  mission  capabilities.  The  main  goal  of  this  thesis  is 
to  provide  a  model,  which  will  be  applied  to  the  current  OUAS  study. 

The  model  allows  the  user  to  define  all  elements  of  the  UAS  via  Excel 
spreadsheet.  The  schedule  generator  then  takes  these  elements  and  generates  a  user 
defined  number  of  individual  schedules  for  every  UAV  included  in  the  scenario.  This 
schedule  generation  is  done  via  Excel  and  Visual  Basic.  Once  the  predefined  number  of 
schedules  is  generated,  the  data  is  read  into  a  General  Algebraic  Modeling  System 
(GAMS)  program  to  find  a  near  optimal  combination  of  individual  schedules. 

B.  AN  INTEGER  PROGRAM  TO  OPTIMIZE  UAS  SCHEDULING 

The  following  integer  linear  program  (ILP),  UAS-ST,  attempts  to  find  the 
consolidated  UAS  schedule  with  the  highest  overall  total  value.  First  the  model  is 
presented,  then  the  detailed  input  required  to  optimize  the  objective  function  is  discussed. 
Once  this  is  complete,  instances  of  UAS-ST  are  generated  to  demonstrate  its  function. 

1.  Indices 


veV 

Tier  levels  [3] 

h^H 

Mission  types  [2]  (H=  {passive,  active }) 

geG 

Ground  Control  Station  (GCS)  [—10] 

1  E  L 

Launch  and  Recovery  Site  (LRS)  [~6] 

m  e  M 

Missions  [-150]  (alias  m’) 
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p  e  P 


Payload  module  types  [—10] 


t  e  T  Time  periods  [96] 

s  e  S  UAV  employment  schedules 

m  ’  €  Pm  Mission  prerequisites:  mission  m  cannot  be  covered  in  any  time 

period  t  unless  at  least  one  mission  m  ’  in  Pm  is  also  covered  in  time 
period  t 

( m,m')eE  Pairwise  exclusive  missions:  mission  m  cannot  be  covered  in  any 
time  period  that  mission  m  ’  is  covered.  E  c:  M  xM  . 

e  A  schedule  s  covers  mission  m  in  time  period  t 

( s,v,t)<=B  schedule  5  uses  UAV  in  tier  level  v  in  time  period  t 

[s,l,t)eK  schedule  5  uses  LRS  /  in  time  period  t 

(, s,c,t )  e  Q  schedule  5  carries  payload  of  class  c  in  time  period  t 

( s,!)<eSL  schedule  5  uses  LRS  / 

(/«, g)  e  MG  mission  m  can  be  covered  by  GCS  g 

( m,h  )  e  MH  mission  m  is  of  mission  type  h 

(g,v)  e  GV  GCS  g  can  support  a  UAS  in  tier  v 

2.  Given  Data  [units] 

valm  Value  per  time  period  of  mission  m  (=total  value  of  m  divided  by 

length  of  m) 

rm  Number  of  time  periods  of  required  coverage  for  mission  m 

num  uavsvj  Number  of  UAVs  in  tier  level  v  available  in  time  period  t 

num _payc,t  Number  of  payloads  of  class  c  available  in  time  period  t 
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gcs  _  cap'l  t  Capacity  (in  UAVs)  of  GCS  g  in  time  period  t  for  missions  of 
type  h 

lrs_capn  Capacity  (in  UAVs)  of  LRS  /  in  time  period  t 
length,,,  Maximum  number  of  periods  mission  m  can  be  covered 

k,  st  =1  if  schedule  5  uses  LRS  /  in  time  period  t 

3.  Decision  Variables 

Xs  =1  if  schedule  5  flown  by  some  UAS  [binary] 

Ymvt  =1  if  mission  m  covered  by  a  tier  v  UAS  in  period  t 

Wmv  =1  if  mission  m  supported  by  GCS  g  for  a  tier  v  UAS  in  time 

period  t  [binary] 

Dm  Total  dwell  time  on  mission  m  [time  periods] 

LOAD.  v ,  Payloads  of  class  c  for  tier  v  sited  at  LRS  I  [cardinality] 
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4. 


Formulation 


max  V  valD„ 

m  m 

m 

(RO) 

s.t.  Z  Xs  <  num _uavslvt 

s:(s,v,t)eB 

(s,l)eSL 

\/l,V,t 

(Rl) 

E  Xs  <  LOADC  VJ 

s:(s,v,t)cB 

(s,c,t)eQ 

\/c,v,l,t 

(R2a) 

YjLOADc,v,i  -  rmm  pay c  v 

l 

Vc,v 

(R26) 

Z  wm,v,g^scs_caphgt 

m,v:(m,g)eMG 

(m,h)eMH 

(g,v)eGK 

\/g,t,h 

(R3) 

X  y  <lrs_caplt 

s:(s,l,t)eK 

VI, t 

(R4) 

y  <  y  x 

m,v,t  —  /  j  s 

s:(s,m,t}eA 

\/m,v,t 

(R5) 

D  <YY  . 

m  m,v,t 

vJ 

Vm 

(R6) 

WliT 

5 

VI 

> 

Vm  :  Pm  ^  0,Vv,Vt 

(R7a) 

y  <  y  w 

m,v,t  m,v,g,t 

g\(m,g)eMG 
(g,v)s  GV 

Vm  :  Pm  =  0,Vv,Vt 

(R76) 

VI 

> 

!>hS 

IXh 

+ 

y 

V(m,m')  e  is,Vt 

(R8) 

D  >r 

m  m 

Vm 

(R9) 

0  <  D.  <  length „ 

Vm 

(RIO) 

X,e{0,l} 

Vs 

(Rl  1) 

Vm,v,t 

(R12) 

Vm,v,g,t 

(R13) 
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5. 


Discussion 


The  objective  function  (RO)  calculates  the  value  of  all  mission  time  periods 
covered.  Constraints  (Rl)  limit  the  number  of  active  UAVs  in  each  tier,  in  each  time 
period,  based  at  each  LRS,  by  the  number  of  UAVS  available.  Constraints  (R2a)  limit  the 
number  of  a  specific  payload  class  an  tier  flown  from  an  LRS  in  each  period  to  the 
number  of  payloads  of  that  class  and  tier  assigned  to  that  LRS,  and  constraints  (R2b) 
limit  the  total  number  of  payloads  of  each  class  and  each  tier  assigned  over  all  LRSs  by 
the  total  number  available  in  that  class.  Constraints  (R3)  limit  the  number  of  missions, 
by  type  supported  by  a  GCS  in  each  time  period.  Constraints  (R4)  limit  the  total  number 
of  UAVs  launching  or  recovering  at  a  LRS  by  the  capacity  of  the  LRS  in  a  given  time 
period.  Constraints  (R5)  controls  whether  or  not  a  mission  is  accomplished  in  any  time 
period.  Constraints  (R6)  limit  the  total  time  spent  on  a  mission  to  the  number  of  time 
periods  covered.  Constraints  (R7a)  address  prerequisite  missions  such  as  AV  CLR 
required  for  long  range  Tier  III  missions,  and  (R7b)  address  line-of-sight  issues  from 
each  GCS  for  missions  that  do  not  have  other  prerequisites.  Constraints  (R8)  prevent 
scheduling  of  mutually  exclusive  missions  by  UAS  of  any  tier.  Constraints  (R9)  force 
“required”  missions  to  be  covered  for  the  number  of  time  periods  required.  Constraints 
(RIO)  limit  the  total  dwell  time  on  a  mission  to  between  zero  and  the  total  amount  of  time 
requested  for  the  mission.  Constraints  (R11-R13)  restrict  schedule  assignment,  mission 
coverage,  and  GCS  assignment  to  be  binary  decisions. 

6.  Route  Enumeration  and  Data  Development 

UAS-ST  requires  a  significant  amount  of  data  processing  to  produce  an  optimal 
schedule;  most  of  this  effort  is  performed  by  a  stack-based  enumeration  routine  that 
generates  all  feasible  schedules  for  each  tier  and  LRS  in  the  given  scenario.  The  schedule 
generator,  or  simply  generator,  requires  two  primary  data  structures  for  computation  of 
feasible  schedules:  a  stack,  PATH(),  containing  a  current  feasible  list  of  missions  to  be 
accomplished,  and  an  array,  ON_STACK(),  that  indicates  whether  each  mission  is 
currently  on  the  stack  of  missions.  Various  data  tables  are  pre-computed  to  aid  in 
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determining  feasibility,  such  as  flight  times  between  airfields  and  missions,  and  between 
pairs  of  missions,  payload  requirements  by  mission,  etc. 

For  each  tier,  v,  and  LRS,  /,  such  that  there  is  at  least  one  UAS  of  tier  v  located  at 
/,  the  generator  starts  with  an  empty  mission  list  PATH(),  (representing  the  trivial  action 
of  launching,  and  then  immediately  recovering,  a  UAS  of  tier  v  from  /),  and  builds 
feasible  mission  lists  exhaustively  by  adding  one  mission  at  a  time  to  the  end  of  the 
current  mission  list  represented  in  PATH().  A  mission  is  only  added  to  the  end  of 
PATH()  if  it  is  within  range  of  a  GCS  capable  of  controlling  the  appropriate  tier,  if  the 
first  period  of  the  mission  that  can  be  accomplished  (based  on  the  current  mission  list) 
allows  the  UAS  to  return  back  to  its  LRS  before  running  out  of  operational  time 
available,  and  if  the  required  list  of  payloads  to  accomplish  all  the  missions  in  PATH() 
can  fit  on  one  UAS  of  tier  v.  If  a  mission  does  not  meet  these  requirements,  it  is 
discarded  and  the  next  mission  from  the  overall  mission  list  is  considered.  Other  rules  for 
calculating  feasible  extensions  to  a  current  feasible  PATH()  can  be  incorporated  easily, 
but  these  three  capture  the  primary  constraints  on  feasibility.  For  example,  we  also  check 
to  see  if  a  mission  is  already  on  the  stack  PATH(),  and  discard  those  to  prevent  missions 
being  revisited.  However,  if  the  capability  to  revisit  missions  in  the  same  schedule  is 
required,  this  test  can  be  removed.  Of  course,  more  feasible  schedules  will  be  generated 
if  this  is  done. 

Once  a  mission  is  added  to  PATH(),  all  missions  are  considered  again  for  the  next 
empty  slot  at  the  end  of  PATH().  In  this  manner  the  generator  provides  a  depth- first 
exploration  of  all  feasible  mission  lists.  If  all  missions  have  been  ruled  out  for  a  slot,  the 
stack  is  “popped,”  the  previous  slot  is  then  considered  again,  and  the  next  mission  in  the 
list  takes  the  current  mission’s  place;  this  replacement  is  repeated  until  a  feasible  mission 
for  that  slot  is  found.  If  no  more  feasible  missions  are  found  for  that  slot,  we  “pop”  the 
stack  again.  When  the  first  slot  is  finished,  we  have  enumerated  all  missions  from  /  using 
tier  v  assets. 

Every  feasible  set  of  missions  in  PATH()  is  recorded  by  incrementing  the  number 
of  feasible  paths  found,  s,  and  then  adding  the  appropriate  elements  to  the  sets 
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[s,m,t)eA,  ( s,v,t)eB ,  ( s,c,t)eQ ,  and  (s,l)eSL  to  define  schedule  s.  See 
Appendix  B  for  the  details  of  the  algorithm. 

We  set  an  upper  bound  on  the  total  number  of  schedules  that  can  be  generated,  so 
as  not  to  create  unsolvable  ILPs.  We  also  provide  a  limit  on  the  number  of  missions 
enumerated  per  schedule,  and  the  generator  implicitly  calculates  a  limit  on  the  number  of 
schedules  to  generate  for  each  tier-LRS  pair,  to  allow  for  the  generation  of  a  nonempty 
list  for  each  such  pair  regardless  of  the  order  in  which  schedules  are  generated  for  those 
pairs.  So,  for  example,  we  might  set  a  limit  of  120,000  total  schedules,  and  four  missions 
per  schedule.  The  generator  calculates  how  many  tier-LRS  pairs  are  feasible,  and  if,  say, 
there  are  six  such  pairs  in  the  scenario,  then  the  generator  will  generate  no  more  than 
20,000  schedules  for  each  pair  considered,  calculated  cumulatively.  (Specifically,  in  this 
example  if  the  first  pair  does  not  generate  all  20,000  of  its  allotted  schedules,  then  the 
generator  can  generate  up  to  40,000  total  schedules  for  the  first  two  pairs  combined.) 
When  these  enumeration  limits  are  reached  for  a  tier-LRS  pair,  the  stack  is  cleared  out 
and  the  generator  moves  to  the  next  tier-LRS  pair. 

This  procedure  can  generate  tens-  or  hundreds-of-thousands  of  feasible  schedules. 
Each  such  schedule  is  associated  with  a  decision  variable,  Xs,  in  the  ILP.  The  limits  on 
enumeration  therefore  restrict  the  schedules  available,  and,  consequently,  the  resulting 
ILP  will  be  a  restriction  of  the  model  that  considers  every  feasible  schedule.  The  more 
schedules  that  are  generated  (through  increasing  the  limit  on  missions  per  schedule,  or  the 
total  schedules  generated,  etc.)  the  better  the  final  schedule  found  by  the  integer  program. 
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IV.  COMPUTATIONAL  RESULTS 


The  OUAS  study  generic  unclassified  scenario  is  replicated  with  the  intent  to  test 
the  results  provided  by  UAS-ST  [OUAS,  2007].  All  computations  are  perfonned  on  a 
Pentium  4  desktop  computer  at  the  Naval  Postgraduate  School  with  the  use  of  the  GAMS 
solver  CPLEX  [2008]. 

A.  INITIAL  TEST  SCENARIO 

The  initial  test  scenario  replicates  a  generic  scenario  given  in  the  OUAS  study 
[OUAS,  2007].  This  replication  provides  the  simplest  and  most  direct  comparison  of 
results  with  ASC-U.  The  scenario  time  period  is  twenty  four  hours  divided  into  ninety 
six  intervals  of  fifteen  minutes  each.  This  scenario  does  not  refer  to  a  specific  country, 
but  is  designed  solely  to  provide  the  framework  for  determining  the  value  of  changes  in 
UAS  structure  on  a  given  set  of  mission  requests.  UAS  force  structure  for  the  initial 
scenario  represents  a  small  but  realistic  composition  of  systems.  Tier  I  is  comprised  of 
three  separate  units  with  four  UAVs  and  one  Ground  Control  Station  (GCS)  per  unit. 
Tier  II  is  comprised  of  three  separate  units  with  one  UAV  and  one  GCS  per  unit.  Tier  III 
is  comprised  of  a  single  unit  with  two  UAVs  and  two  GCSs. 

Table  2  provides  a  summary  of  the  given  UAS  force  structure  for  the  generic  scenario. 


Number  of: 

Tier  I 

Tier  II 

Tier  III 

Systems 

12 

3 

2 

UAVs  per  system 

3 

3 

4 

GCSs  per  system 

1 

1 

2 

UAVs  controlled  by  GCS 

1 

2 

2 

Table  2.  Number  of  UAS  Assets  Available  in  Generic  Scenario 


The  initial  test  scenario  consists  of  placing  units  from  different  UAS  Tiers  at 
specific  ranges  to  test  the  scheduling  constraints  of  the  model.  The  tier  III  squadron  is 

centered  at  an  airfield  in  support  of  a  regiment.  The  squadron  consists  of  two  UAVs  and 
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two  GCSs  at  this  location.  These  are  the  only  tier  III  assets  available  for  our  test 
scenario.  The  mission  areas  for  the  tier  III  UAS  are  obtained  by  plotting  points  in  the 
cardinal  directions  at  or  beyond  one  hundred  and  fifty  nautical  miles,  (beyond  tier  II 
range). 

The  three  tier  II  units  are  deployed  in  a  triangular  pattern  around  the  tier  III 
location  in  direct  support  of  individual  battalions.  The  tier  II  units  are  separated  from 
each  other  by  a  distance  of  at  least  20  nautical  miles.  As  before,  specific  mission  areas 
are  chosen  for  each  tier  II  unit. 

Tier  I  units  are  collocated  in  the  same  manner  as  tier  II,  but  are  in  direct  support  at 
the  company  level.  Each  tier  I  unit  consists  of  four  UAVs  and  a  single  GCS.  Mission 
areas  are  selected  with  a  special  emphasis  on  their  range.  Once  the  areas  are  designated, 
they  are  used  repetitively  to  generate  multiple  missions.  The  use  of  repetitive  mission 
areas  makes  it  much  simpler  to  detect  mission  assignments  beyond  the  range  of  a  specific 
tier.  Tier  I  and  II  missions  are  within  range  of  tier  III  assets. 


lTash,l 


Figure  3.  Plot  of  OUAS  Scenario  Taken  from  OUAS 
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B. 


UAS-ST  INITIAL  RESULTS  AND  ANALYSIS 


Data  for  the  initial  scenario  is  entered  via  a  mission  request  worksheet  in  Excel. 
The  user  enters  mission  requests  for  UAS  coverage  in  the  following  manner: 


Mission  ID 

Location 

Mission 

Class 

Start 

Time 

End 

Time 

Total 

Time 

Requested 

Start 

Period 

Finish 

Period 

Time 

Periods 

Requested 

Value 

Per  Time 

Period 

Lat 

Long 

ml 

N323934 

W1 144850 

AV  CLR 

0:00 

8:00 

8:00 

1 

32 

32 

5 

m2 

N322947 

W1 144606 

ISR 

0:30 

2:30 

2:00 

3 

10 

8 

10 

m3 

N322825 

W1 144547 

ISR 

1:00 

4:00 

3:00 

5 

16 

12 

10 

m4 

N322648 

W1 144730 

ISR 

1:30 

2:30 

1:00 

7 

10 

4 

10 

m5 

N322549 

W1 144914 

ISR 

2:00 

4:00 

2:00 

9 

16 

8 

10 

m6 

N324012 

W1 142426 

ISR 

2:30 

5:30 

3:00 

11 

22 

12 

10 

m7 

N324021 

W1 142011 

ISR 

3:00 

4:00 

1:00 

13 

16 

4 

10 

m8 

N324055 

W1 142318 

ISR 

3:30 

5:30 

2:00 

15 

22 

8 

10 

m9 

N323922 

W1 142442 

ISR 

4:00 

7:00 

3:00 

17 

28 

12 

10 

mIO 

N323848 

W1 142056 

ISR 

4:30 

5:30 

1:00 

19 

22 

4 

10 

Table  3.  Input  Worksheet  for  Mission  Requests 


Waypoints  for  the  mission  areas  are  placed  at  specific  ranges  to  test  the 
scheduling  constraints  of  the  model.  These  waypoints  are  then  used  to  generate  a  table  of 
one  hundred  and  fifty  mission  requests.  The  requests  are  broken  down  by  mission  areas. 
Mission  requests  one  through  forty  focus  on  the  tier  I  mission  areas.  Requests  forty  one 
through  one  hundred  and  three  focuses  on  tier  II  and  the  remaining  requests  are  tier  III 
mission  areas.  See  Appendix  A  for  the  complete  scenario  worksheet  used  in  analysis  of 
UAS-ST.  The  model  constraints  allow  a  higher  tier  vehicle  to  perfonn  a  lower  tier 
mission  if  that  is  the  highest  value  mission  available  at  the  time  and  the  vehicle  is  within 
range.  However,  the  longer  range  of  each  subsequent  tier’s  mission  area  make  it 
infeasible  to  schedule  a  lower  tier  UAS.  The  sequential  mission  numbers  enable  the  user 
to  quickly  identify  any  infeasible  mission  assignments. 

Mission  classes  are  defined  according  to  the  mission  categories  given  in  the 

OUAS  study.  Values  for  the  individual  missions  are  chosen  arbitrarily.  Selection  of  the 

assigned  mission  values  is  critical  because  they  have  the  ability  to  skew  the  results  by 

inflating  the  value  of  a  specific  category.  In  the  scenario,  active  missions  are  given 

greater  priority  than  passive.  Of  the  active  subset,  missions  which  are  unique  to  a 
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specific  tier  are  given  the  highest  priority.  Missions  which  are  redundant  to  all  tiers  are 
given  the  lowest  priority.  Table  4  outlines  the  missions  which  each  tier  can  perform  and 
their  associated  value  in  the  scenario. 


Tier  1 

Tier  II 

Tier  III 

ISR 

10 

15 

20 

AVCLR 

5 

10 

5 

CR 

10 

5 

LSR  P 

30 

SIGINT 

30 

EW 

40 

STK 

60 

LSR  D 

40 

Table  4.  Assigned  Mission  Values  by  Tier 


The  requested  mission  time  is  defined  through  the  start  and  end  times,  which  are 
entered  by  the  user.  UAS-ST  then  converts  that  time  into  the  appropriate  number  of 
fifteen  minute  time  periods.  Once  a  UAS  is  assigned  in  a  given  time  period  it  may  not  be 
reassigned  until  the  next  time  period.  Times  are  started  at  the  beginning  of  the  first  time 
period  and  continue  to  the  end  of  the  final  period. 

Another  element  which  greatly  affects  the  optimal  value  achieved  is  the  number 
of  mission  payloads  available.  See  Figure  4  for  a  sample  of  the  payload  worksheet.  By 
manipulating  the  number  of  payloads  of  a  given  type  available  to  a  UAS  tier,  it  is 
possible  to  test  the  scheduling  constraints  and  see  very  quickly  if  payload  availability  was 
a  limiting  factor.  This  is  highly  beneficial  in  detennining  future  requirements  for  UAS 
structure. 
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Payload 

ID 

Payload  Type 

Mission 

Class 

Quantity 

Available 

UAS 

Tier 

Pi 

EO/IR 

ISR 

9 

T1 

P2 

EO/IR-targeting 

ISR  T 

3 

T2 

p3 

EO/IR-targeting 

ISR  T 

2 

T3 

P4 

SAR  w/MTI 

ISR  SAR 

1 

T3 

p5 

SIGINT 

SIGINT 

1 

T2 

P6 

SIGINT 

SIGINT 

1 

T3 

P7 

Mine  Detection 

ISR  MD 

1 

T3 

P8 

Pointer 

LSR  P 

1 

T2 

p9 

Pointer 

LSR  P 

1 

T3 

plO 

Rangefinder 

LSR  RF 

1 

T3 

pi  1 

Designator 

LSR  D 

1 

T3 

p12 

AV  CLR 

AV  CLR 

9 

T1 

p13 

AV  CLR 

AV  CLR 

3 

T2 

p14 

AV  CLR 

AV  CLR 

1 

T3 

p15 

CR 

CR 

1 

T2 

p16 

CR 

CR 

1 

T3 

p17 

Strike 

STK 

2 

T3 

p18 

EW 

EW 

2 

T3 

Figure  4.  Example  of  Payload  Worksheet 


System  capabilities  and  limitations  are  drawn  directly  from  the  OUAS  study.  The 
study  seeks  to  define  UAS  tier  requirements.  Therefore,  those  assumptions  are  critical  to 
any  analysis  that  is  conducted.  The  tier  parameters  as  given  in  Figure  5  are  used  in  all 
testing.  Operational  time  limit  and  maintenance  time  are  not  known  at  this  time.  Those 
columns  are  for  future  expansion  when  the  data  becomes  available. 


UAS  Tier 

Max 

Endurance  (hrs) 

Cruise 

Speed(knots) 

Operational 

Radius(nm) 

Payload 

Capacity 

Operational 
Time  Limit(hrs) 

Maintenance 

Time 

T 1 

1.5 

17 

5 

1 

T2 

15 

60 

50 

2 

T3 

8 

200 

300 

3 

*  All  Values  Taken  from  Table  ES-2  of  Overarching  UAS  Study  21Nov2007 

Figure  5.  OUAS  Study  Future  Tier  Performance  Parameters 

Once  all  data  for  the  scenario  is  entered  via  Excel  spreadsheet,  UAS-ST  is  run 
from  the  graphical  interface  given  in  Figure  6.  The  interface,  known  as  the  dashboard, 
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allows  the  user  to  specify  several  key  parameters.  First  the  user  specifies  the  start  time, 
the  time  periods  in  horizon,  and  the  number  of  periods  per  hour.  For  the  mission 
generation  portion,  a  user  specifies  the  maximum  number  of  schedules  that  UAS-ST  can 
generate  and  the  maximum  number  of  missions  per  schedule  that  a  single  UAV  can 
perform.  These  two  parameters  greatly  affect  the  outcome  of  the  model. 

As  the  number  of  maximum  missions  per  schedule  increases,  the  number  of 
possible  schedules  rapidly  increases  as  well.  Setting  the  number  of  maximum  schedules 
too  low  prevents  a  full  enumeration  of  all  possible  schedules  and  is  a  constraint  on  the 
model.  Setting  the  number  of  maximum  schedules  too  high  allows  a  full  enumeration  of 
all  possible  schedules  and  results  in  a  longer  processing  time.  See  Figure  6  for  an 
example.  All  tested  combinations  satisfy  the  requirement  for  a  quick  total  solution  time. 
Once  all  parameters  are  set,  selecting  the  BUILD  function  initiates  VBA  to  build  the 
designated  number  of  flight  schedules.  After  the  build  is  complete,  selecting  SOLVE 
initiates  GAMS  to  read  in  the  data,  apply  the  selected  solver,  and  generate  a  near  optimal 
solution.  When  GAMS  reaches  a  solution,  selection  of  RESULTS  displays  the  solution  in 
the  format  seen  in  Figure  7. 


UAS  Planner  vl.01 

Start  Time 

0:00 

Time  Periods  in  Horizon 

96 

First  Planning  Period 

1 

Last  Planning  Period 

96 

Periods  per  hour 

4 

Max  Schedules 

75000 

Max  Missions  per  Schedule 

3 

Solver 

CPLEX 

OPTCR 

0.05 

RESLIM 

3600 

2-Jun-08 


Build 


Solve 


Results 


Routes: 

Status: 


67324 


Figure  6.  Sample  of  User  Interface  (Dashboard) 
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Total 

Total 

Value 

Total 

Start 

End 

Time  Periods 

Time  Periods 

Per  Time 

Value 

Mission 

Time 

Time 

Requested 

Assigned 

Period 

Achieved 

ml 

0:00 

8:00 

32 

0 

10 

0 

m2 

0:30 

2:30 

8 

0 

5 

0 

m3 

1:00 

4:00 

12 

0 

5 

0 

m4 

1:30 

2:30 

4 

0 

5 

0 

m5 

2:00 

4:00 

8 

0 

5 

0 

m6 

2:30 

5:30 

12 

6 

5 

30 

ml 

3:00 

4:00 

4 

1 

5 

5 

m8 

3:30 

5:30 

8 

4 

5 

20 

m9 

4:00 

7:00 

12 

7 

5 

35 

mIO 

4:30 

5:30 

4 

0 

5 

0 

ml  1 

5:00 

7:00 

8 

5 

5 

25 

m12 

5:30 

8:30 

12 

9 

5 

45 

m13 

6:00 

7:00 

4 

0 

5 

0 

m14 

6:30 

8:30 

8 

0 

5 

0 

m15 

7:00 

10:00 

12 

6 

5 

30 

m16 

7:30 

8:30 

4 

0 

5 

0 

m17 

8:00 

16:00 

32 

0 

10 
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Figure  7.  Partial  Display  of  Optimal  UAS-ST  Schedule 


Figure  7  demonstrates  the  format  in  which  the  solution  given  by  UAS-ST  is 
displayed.  All  missions  are  sorted  sequentially  by  their  assigned  mission  number.  The 
total  value  achieved  for  each  mission  is  given  in  the  right  column.  The  total  value 
obtained  by  the  entire  schedule  is  given  in  the  top  right  corner.  This  format  enables  the 
easiest  comparison  between  runs  because  it  quickly  lets  the  user  see  which  missions  are 
scheduled  and  for  what  portion  of  the  requested  time.  If  certain  missions  are  not  being 
scheduled,  the  trend  is  easy  to  detect. 

C.  COMPARING  UAS-ST  TO  ASC-U 

ASC-U  is  the  model  originally  used  in  the  OUAS  study.  It  is  highly  useful  in  the 
analysis  of  UAS  future  requirements  with  several  strong  attributes.  However,  there  are 
some  areas  which  require  modification  to  meet  the  needs  of  the  OUAS  study.  The  intent 
of  this  thesis  is  to  develop  a  user  friendly  model  that  addresses  those  shortcomings  and 
meets  the  requirements  of  OAD. 
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For  ease-of-use  the  individual  elements  of  the  UAS  are  separated  into  different 
worksheets.  A  separate  spreadsheet  for  the  UAS  tier-  specific  performance  requirements 
is  also  included.  This  spreadsheet  allows  the  user  to  verify  the  basic  parameters  of  all 
UAS  tiers  in  one  chart.  The  user  can  easily  modify  a  parameter  and  re-run  the  model  to 
analyze  the  effect.  The  model  applies  these  values  to  the  system  constraints  when 
optimizing  the  schedule. 

UAS-ST  takes  a  different  approach  from  ASC-U  when  addressing  the 
optimization  of  payloads.  In  ASC-U  all  possible  combination  of  payloads  are 
enumerated  and  then  designated  as  mission  packages.  UAS-ST  treats  the  payloads  as 
individual  units,  which  allows  for  tracking  of  specific  payload  utilization.  OAD  indicated 
this  was  necessary  for  determining  future  requirements.  The  payload  data  is  further 
separated  by  creating  a  mission  class  worksheet  and  a  payload  worksheet.  The  mission 
class  worksheet  defines  all  current  missions  for  the  UAS.  It  allows  the  user  to  designate 
a  mission  as  either  active  or  passive  and  which  tier  is  capable  of  covering  that  mission. 
The  payload  worksheet  allows  the  user  to  designate  the  quantity  available. 

A  major  difference  between  the  two  models  is  in  the  handling  of  the  mission  data. 
In  UAS-ST  all  information  is  consolidated,  both  input  and  output,  in  the  same  worksheet. 
The  user  builds  the  missions  page  in  the  same  manner  as  a  flight  schedule  and  assigns  a 
specific  number  to  an  individual  mission.  UAS-ST  sorts  all  output  based  on  mission 
number,  which  makes  for  easy  tracking  of  a  specific  mission.  ASC-U  uses  either  the 
mission  number  or  the  coordinates  of  the  mission  location  to  sort  output.  This  makes 
analysis  of  results  extremely  time  consuming  and  requires  additional  sorting.  ASC-U 
requires  the  use  of  a  heuristic,  SUPER  MISSIONS,  to  designate  a  mission  as  mandatory 
by  assigning  it  extra  value.  UAS-ST  addresses  this  in  the  mission  input  page  by  allowing 
the  user  to  designate  any  portion  of  a  mission  request  as  required.  Setting  this 
requirement  serves  as  an  additional  constraint  on  the  optimization  model. 

ASC-U  utilizes  an  Excel  spreadsheet  and  pulls  the  data  into  an  Access  database. 
UAS-ST  applies  an  Excel  spreadsheet,  which  uses  Visual  Basic  to  transfer  the  data  to 
GAMS.  GAMS  then  applies  the  CPLEX  solver  to  quickly  optimize  the  schedule  for  the 

mission  requests  it  is  given.  The  required  features  addressed  by  the  supplemental 
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heuristics  in  ASC-U  are  incorporated  directly  into  UAS-ST  to  eliminate  the  need  for 
further  processing.  This  allows  for  a  much  quicker  and  efficient  analysis  of  alternatives. 
CPLEX  is  highly  efficient  in  preprocessing  the  feasible  solutions  and  reduces  the  run 
time  to  a  matter  of  minutes.  Overall,  UAS-ST  provides  an  efficient  update  to  ASC-U. 
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V.  CONCLUSIONS  AND  NEW  OPPORTUNITIES 


A.  SUMMARY 

This  thesis  develops  a  scheduling  tool,  UAS-ST  as  a  replacement  for  ASC-U. 
UAS-ST  addresses  several  areas  which  require  additional  processing  in  ASC-U.  UAS- 
ST  allows  a  user  to  generate  a  daily  flight  based  on  a  given  set  of  operational  and 
performance  parameters.  The  true  strength  of  UAS-ST  is  that  a  user  can  change  a  single 
parameter  or  some  combination  of  parameters  and  quickly  rerun  the  model  to  see  the 
effect  on  a  flight  schedule.  UASs  are  critical  to  our  nation’s  future  combat  effectiveness. 
As  a  result,  all  branches  of  DoD  recognize  the  need  to  develop  an  integrated  network  of 
UASs.  MCCDC  OAD  is  currently  conducting  an  OUAS  study  with  the  goal  of  defining 
future  UAS  requirements.  This  thesis  is  a  direct  contribution  to  that  study,  and  the  model 
developed  herein,  provides  a  planning  tool  for  the  development  of  future  UAS  structure. 

B.  OPERATIONAL  INTRODUCTION 

UAS-ST  is  currently  intended  for  use  in  the  second  phase  of  the  OUAS  study.  It 
will  replace  ASC-U  and  provide  another  tool  for  quantitative  analysis.  UAS-ST  provides 
new  capabilities  which  ASC-U  was  unable  to  perfonn.  With  the  use  of  the  CPLEX 
solver  in  GAMS,  UAS-ST  is  able  to  quickly  analyze  tens  of  thousands  of  schedules  for 
individual  UAVs  and  then  choose  a  near-optimal  subset  to  provide  a  complete  schedule. 
The  user  is  able  to  modify  all  aspects  of  the  UAS  structure  through  a  simple  spreadsheet 
interface.  Output  is  provided  in  the  same  format,  which  makes  understanding  the  results 
and  error  tracking  much  simpler.  ASC-U  requires  the  use  of  several  heuristics  to  meet 
the  needs  of  OUAS.  UAS-ST  directly  addresses  these  issues  without  the  need  for 
additional  heuristics.  UAS-ST  provides  a  rapid  analysis  tool  and  should  be  implemented 
immediately  in  the  ongoing  OUAS  study. 
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c. 


FUTURE  DEVELOPMENT 


UAS  development  is  expanding  rapidly  with  several  new  technologies  becoming 
available.  Therefore,  modeling  requirements  are  expected  to  change  rapidly  as  well.  A 
possible  issue  with  UAS-ST  is  the  requirement  for  the  GAMS  CPLEX  solver  to  produce 
a  quick  solution.  Due  to  portability  issues,  the  future  development  of  a  non-proprietary 
heuristic  is  a  valuable  topic  for  future  analysis.  To  improve  on  UAS-ST,  the  heuristic 
needs  to  modify  the  enumeration  routine  used  in  the  model.  Rather  than  a  total 
enumeration  of  all  possible  schedules,  a  new  heuristic  should  initially  select  higher  value 
routes  and  not  consider  low  value  routes  beneath  a  designated  limit.  One  possible 
method  is  a  greedy  heuristic  which  selects  the  highest  value  subset  of  individual 
schedules  and  eliminate  all  others  from  consideration.  By  initially  eliminating  low  value 
routes,  the  overall  processing  time  is  greatly  reduced.  In  the  long  term,  analysts  in  theater 
may  run  UAS-ST  on  their  laptops  as  they  plan  for  actual  UAS  employment.  The  ultimate 
goal  is  to  find  a  near-optimal  solution  without  the  need  for  specialty  software 
applications. 
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APPENDIX  A  TEST  SCENARIO 


Tier  I  Mission  Areas  (ml-m40) 
Tier  II  Mission  Areas  (m41-ml03) 
Tier  III  Mission  Areas  (ml 04- 150) 


Mission  ID 

Location 

Mission 

Class 

Start 

Period 

Finish 

Period 

Time 

Periods 

Requested 

Value 

Per  Time 

Period 

Lat 

Long 

ml 

N323934 

W1 144850 

AV  CLR 

1 

32 

32 

5 

m2 

N322947 

W1 144606 

ISR 

3 

10 

8 

10 

m3 

N322825 

W1 144547 

ISR 

5 

16 

12 

10 

m4 

N322648 

W1 144730 

ISR 

7 

10 

4 

10 

m5 

N322549 

W1 144914 

ISR 

9 

16 

8 

10 

m6 

N324012 

W1 142426 

ISR 

11 

22 

12 

10 

m7 

N324021 

W1 142011 

ISR 

13 

16 

4 

10 

m8 

N324055 

W1 142318 

ISR 

15 

22 

8 

10 

m9 

N323922 

W1 142442 

ISR 

17 

28 

12 

10 

mIO 

N323848 

W1 142056 

ISR 

19 

22 

4 

10 

ml  1 

N325910 

W1 142608 

ISR 

21 

28 

8 

10 

m12 

N325733 

W1 142738 

ISR 

23 

34 

12 

10 

m13 

N325757 

W1 142944 

ISR 

25 

28 

4 

10 

m14 

N325927 

W1 142938 

ISR 

27 

34 

8 

10 

m15 

N330047 

W1 142630 

ISR 

29 

40 

12 

10 

m16 

N323934 

W1 144850 

ISR 

31 

34 

4 

10 

m17 

N322947 

W1 144606 

AV  CLR 

33 

64 

32 

5 

m18 

N322825 

W1 144547 

ISR 

35 

46 

12 

10 

m19 

N322648 

W1 144730 

ISR 

37 

40 

4 

10 

m20 

N322549 

W1 144914 

ISR 

39 

46 

8 

10 

m21 

N324012 

W1 142426 

ISR 

41 

52 

12 

10 

m22 

N324021 

W1 142011 

ISR 

43 

46 

4 

10 

m23 

N324055 

W1 142318 

ISR 

45 

52 

8 

10 

m24 

N323922 

W1 142442 

ISR 

47 

58 

12 

10 

m25 

N323848 

W1 142056 

ISR 

49 

52 

4 

10 

m26 

N325910 

W1 142608 

ISR 

51 

58 

8 

10 

m27 

N325733 

W1 142738 

ISR 

53 

64 

12 

10 

m28 

N325757 

W1 142944 

ISR 

55 

58 

4 

10 

m29 

N325927 

W1 142938 

ISR 

57 

64 

8 

10 

m30 

N330047 

W1 142630 

ISR 

59 

70 

12 

10 

m31 

N323934 

W1 144850 

ISR 

61 

64 

4 

10 

m32 

N324012 

W1 142426 

ISR 

63 

70 

8 

10 

m33 

N325910 

W1 142608 

AV  CLR 

65 

96 

32 

5 

m34 

N322947 

W1 144606 

ISR 

67 

70 

4 

10 

37 


m35 

N324021 

W1 142011 

ISR 

69 

76 

8 

10 

m36 

N325733 

W1 142738 

ISR 

71 

82 

12 

10 

m37 

N322825 

W1 144547 

ISR 

73 

76 

4 

10 

m38 

N324055 

W1 142318 

ISR 

75 

82 

8 

10 

m39 

N325757 

W1 142944 

ISR 

77 

88 

12 

10 

m40 

N322648 

W1 144730 

ISR 

79 

82 

4 

10 

m41 

N321813 

W1 145516 

ISR 

1 

4 

4 

15 

m42 

N321817 

W1 143321 

ISR 

3 

18 

16 

15 

m43 

N322752 

W1 150128 

ISR 

5 

36 

32 

15 

m44 

N322413 

W1 140639 

AV  CLR 

1 

48 

48 

10 

m45 

N324657 

W1 135856 

CR 

1 

48 

48 

10 

m46 

N325149 

W1 142341 

LSR  P 

1 

8 

8 

30 

m47 

N330310 

W1 144304 

ISR 

7 

10 

4 

15 

m48 

N331121 

W1 141606 

ISR 

9 

24 

16 

15 

m49 
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W1 144635 

ISR 

11 

42 

32 

15 
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9 

56 
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10 
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N321817 

W1 143321 

CR 

9 
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10 
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8 

30 
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15 
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48 
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22 
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15 
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30 

m65 

N330310 

W1 144304 

ISR 

25 

28 
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42 
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33 
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64 

8 

30 

m71 

N322413 
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34 
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48 

16 
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35 
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88 
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37 

40 

4 

15 
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W1 143321 
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39 

54 

16 

15 
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41 

72 

32 

15 
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W1 140639 
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49 

96 

48 

10 
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N324657 

W1 135856 
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49 

96 

48 

10 

38 


m82 

N325149 

W1 142341 

LSR  P 

81 

88 

8 

30 

m83 

N330310 

W1 144304 

ISR 

43 

46 

4 

15 

m84 

N331121 

W1 141606 

ISR 

45 

60 

16 

15 

m85 

N331641 

W1 144635 

ISR 

47 

78 

32 

15 

m86 

N321813 

W1 145516 

AV  CLR 

57 

96 

40 

10 

m87 

N321817 

W1 143321 

CR 

57 

96 

40 

10 

m88 

N322752 

W1 150128 

LSR  P 

89 

96 

8 

30 

m89 

N322413 

W1 140639 

ISR 

49 

52 

4 

15 

m90 

N324657 

W1 135856 

ISR 

51 

66 

16 

15 

m91 

N325149 

W1 142341 

ISR 

53 

84 

32 

15 

m92 

N330310 

W1 144304 

ISR 

55 

58 

4 

15 

m93 

N331121 

W1 141606 

ISR 

57 

72 

16 

15 

m94 

N331641 

W1 144635 

ISR 

59 

82 

24 

15 

m95 

N321813 

W1 145516 

ISR 

61 

64 

4 

15 

m96 

N321817 

W1 143321 

ISR 

63 

66 

4 

15 

m97 

N322752 

W1 150128 

ISR 

65 

80 

16 

15 

m98 

N322413 

W1 140639 

ISR 

67 

96 

30 

15 

m99 

N324657 

W1 135856 

ISR 

69 

72 

4 

15 

mlOO 

N325149 

W1 142341 

ISR 

71 

86 

16 

15 

mIOl 

N330310 

W1 144304 

ISR 

73 

96 

24 

15 

ml  02 

N331121 

W1 141606 

ISR 

75 

78 

4 

15 

m103 

N331641 

W1 144635 

ISR 

77 

92 

16 

15 

m104 

N312538 

W1 145510 

ISR 

1 

96 

96 

20 

m105 

N334259 

W1 155906 

ISR 

1 

48 

48 

20 

m106 

N333540 

Will  5204 
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49 

96 

48 

20 

m107 

N312538 

W1 145510 
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1 

32 

32 

20 

m108 

N334259 
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33 

64 

32 

20 

m109 

N333540 
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65 

96 

32 

20 

mllO 

N312538 

W1 145510 

CR 

1 

96 

96 

5 

mil  1 

N334259 

W1 155906 

CR 

1 

96 

96 

5 

ml  12 

N333540 

Will  5204 

CR 

1 

96 

96 

5 

ml  13 

N312538 

W1 145510 

CR 

1 

96 

96 

5 

ml  14 

N334259 

W1 155906 

SIGINT 

1 

96 

96 

30 

ml  15 

N333540 

W1 115204 

SIGINT 

1 

96 

96 

30 

ml  16 

N312538 

W1 145510 

SIGINT 

1 

96 

96 

30 

ml  17 

N334259 

W1 155906 

SIGINT 

1 

96 

96 

30 

m118 

N333540 

Will  5204 

EW 

1 

8 

8 

40 

m119 

N312538 

W1 145510 

EW 

9 

16 

8 

40 

m120 

N334259 

W1 155906 

EW 

25 

32 

8 

40 

m121 

N333540 

Will  5204 

EW 

33 

40 

8 

40 

ml  22 

N312538 

W1 145510 

EW 

49 

56 

8 

40 

m123 

N334259 

W1 155906 

EW 

57 

64 

8 

40 

ml  24 

N333540 
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EW 

73 

80 

8 

40 
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W1 145510 

EW 

81 

88 

8 

40 
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W1 155906 
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1 

16 

16 

60 

m127 
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Will  5204 

STK 

17 

32 

16 

60 

m128 

N312538 

W1 145510 

STK 

33 

48 

16 

60 

39 


m129 

N334259 

W1 155906 

STK 

49 

64 

16 

60 

m130 

N333540 

Will  5204 

STK 

65 

80 

16 

60 

m131 

N312538 

W1 145510 

STK 

81 

96 

16 

60 

ml  32 

N334259 

W1 155906 

STK 

25 

40 

16 

60 

m133 

N333540 

Will  5204 

STK 

73 

88 

16 

60 

ml  34 

N312538 

W1 145510 

LSR  D 

1 

8 

8 

40 

m135 

N334259 

W1 155906 

LSR  D 

25 

32 

8 

40 

m136 

N333540 

Will  5204 

LSR  D 

49 

56 

8 

40 

m137 
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73 

80 

8 

40 
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24 

24 

5 
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48 

24 

5 
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72 

24 

5 
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73 

96 

24 

5 
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96 
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5 
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24 
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APPENDIX  B  VBA  CODE 


Sub  SolveUAS ( ) 

Dim  bSolveResult  As  Boolean 
bSolveResult  =  SolveProblem (True) 

End  Sub 

Sub  BuildData ( ) 

Dim  uasDataFN  As  Integer 
Dim  uasSetsFN  As  Integer 
Dim  uasDynSetsFN  As  Integer 
Dim  logFN  As  Integer 
Dim  rTiers  As  Range 
Dim  rGCSs  As  Range 
Dim  rLRSs  As  Range 
Dim  rUAVs  As  Range 
Dim  rMissions  As  Range 
Dim  rPayloads  As  Range 
Dim  rSchedules  As  Range 
Dim  rClasses  As  Range 

Dim  INumTiers  As  Long,  INumGCSs  As  Long,  INumLRSs  As  Long,  INumUAVs  As 

Long 

Dim  INumMissions  As  Long,  INumPayloads  As  Long,  INumClasses  As  Long 
Dim  INumSchedules  As  Long 
Dim  INumLRSSchedules  As  Long 
Dim  bNextRoute  As  Boolean 

Dim  rvalues  As  Range 
Dim  INumValues  As  Long 
Dim  rDependencies  As  Range 
Dim  lPlanPeriods  As  Long 
Dim  lFirstPeriod  As  Long 
Dim  lLastPeriod  As  Long 
Dim  IPeriodsPerHour  As  Long 

Dim  ICount  As  Long 

Dim  currCell  As  Range 

Dim  iRow  As  Integer,  iCol  As  Integer 

Dim  IMission  As  Long,  !Mission2  As  Long,  1UAV  As  Long 
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Dim  IPeriod  As  Long,  lPeriodl  As  Long,  lPeriod2  As  Long 

Dim  ITier  As  Long,  IPayload  As  Long,  1GCS  As  Long,  1LRS  As  Long 

Dim  lClass  As  Long 

Dim  tierNames  As  StringDictionaryClass 
Dim  uavNames  As  StringDictionaryClass 
Dim  gcsNames  As  StringDictionaryClass 
Dim  IrsNames  As  StringDictionaryClass 
Dim  missionNames  As  StringDictionaryClass 
Dim  payloadNames  As  StringDictionaryClass 
Dim  classNames  As  StringDictionaryClass 
Dim  periodNames  As  StringDictionaryClass 
Dim  tValues  As  TupleDictionaryClass 


uasDataFN  =  FreeFileO 
Open  ThisWorkbook . path  & 
uasSetsFN  =  FreeFileO 
Open  ThisWorkbook . path  & 
logFN  =  FreeFileO 
Open  ThisWorkbook . path  & 


"\uasData . gms"  For  Output  As  uasDataFN 


"\uasSets . gms"  For  Output  As  uasSetsFN 


"\uasXL.log"  For  Output  As  logFN 


{STATIC,  GROUNDED}  SETS 


v 

h 

g 

1 


assigned} 


t 

P 


s 


Tier  levels 

mission  types  {active,  passive} 

Ground  control  stations 
Launch  and  recovery  sites 

Missions  {missions  that  ships  can  be 

Time  periods 
Payload  module  types 
UAV  employment  schedules 


Print 

#logFN,  "Starting  output" 

Print 

#uasSetsFN, 

II 

OPTIONS" 

Print 

#uasSetsFN, 

II 

MIP  =  "  & 

Print 

#uasSetsFN, 

II 

OPTCR  =  " 

Print 

#uasSetsFN, 

II 

RESLIM  = 

Print 

#uasSetsFN, 

II 

.  II 

wsDashboard . Range ("SOLVER") 

&  wsDashboard . Range ( "OPTCR" ) 

"  &  wsDashboard. Range ( "RESLIM" ) 
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Print  #uasSetsFN,  "  SETS" 

Set  rTiers  =  wsTiers . Range ( "A4 " ,  wsTiers . Range ( "A4" ). End (xlDown) ) 

Call  GAMSWr iteSetDef ( rTiers ,  "v",  uasSetsFN,  "Tier  levels") 

Set  rUAVs  =  wsUAV . Range ( "A4 " ,  wsUAV . Range ( "A4 "). End (xlDown) ) 

Print  #uasSetsFN,  "  h  /" 

Print  #uasSetsFN,  "  h_active" 

Print  #uasSetsFN,  "  h_passive" 

Print  #uasSetsFN,  "  /" 

Set  rGCSs  =  wsGCS . Range ( "A4 " ,  wsGCS . Range ( "A4 "). End (xlDown) ) 

Call  GAMSWr iteSetDef (rGCSs,  "g",  uasSetsFN,  "Ground  control  stations") 
Set  rLRSs  =  wsLRS . Range ( "A4 " ,  wsLRS . Range ( "A4 "). End (xlDown) ) 

Call  GAMSWriteSetDef (rLRSs,  "1",  uasSetsFN,  "Launch  and  recovery 

sites") 

Set  rMissions  =  wsMissions . Range ( "A4 " , 

wsMissions . Range ( "A4 " ) .End (xlDown) ) 

Call  GAMSWriteSetDef (rMissions,  "m" ,  uasSetsFN,  "Missions") 

Set  rPayloads  =  wsPayloads . Range ( "A4 " , 

ws Pay loads . Range ( "A4 " ) .End (xlDown) ) 

'  Call  GAMSWriteSetDef (rPayloads ,  "p",  uasSetsFN,  "Payloads") 

Set  rClasses  =  wsMClasses . Range ( "A4 " , 

wsMClasses . Range ( "A4" ) . End (xlDown) ) 

Call  GAMSWriteSetDef (rClasses,  "c",  uasSetsFN,  "Payload  Classes") 

IPeriodsPerHour  =  CLng (wsDashboard. Range ( "PERIODS_PER_HOUR" ) ) 
lPlanPeriods  =  CLng  (wsDashboard . Range ( "HORI ZON" ) ) 
lFirstPeriod  =  CLng (wsDashboard . Range ( "FIRST_PERIOD" ) ) 
lLastPeriod  =  CLng (wsDashboard. Range ( "LAST_PERIOD" ) ) 

If  lLastPeriod  >  lPlanPeriods  Then 

MsgBox  "Last  Period  comes  after  maximum  planning  horizon: 

truncating . " 

lLastPeriod  =  lPlanPeriods 
End  If 

If  lFirstPeriod  >  lLastPeriod  Then 

MsgBox  "Empty  planning  horizon:  first  Period  occurs  after  last 
Period.  Aborting." 

GoTo  subExit 

End  If 

Print  #logFN,  "Periods  "  &  lPlanPeriods  &  "  First:  "  &  lFirstPeriod  & 
"  Last:  "  &  lLastPeriod 

Print  fuasSetsFN,  "  t  Time  Periods  /t_"  &  Format  ( lFirstPeriod)  & 

"*t  "  &  Format ( lLastPeriod)  &  "/" 
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Set  tierNames  =  New  Str ingDictionaryClass 

If  Def ineGroundedSet ( rTiers ,  tierNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (tiers) .  User  chose 
to  cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumTiers  =  tierNames . Count 

For  ICount  =  1  To  INumTiers 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  tierNames . Item (ICount) 

Next  ICount 

Set  uavNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet (rUAVs,  uavNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (uavs)  .  User  chose 
to  cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumUAVs  =  uavNames . Count 

For  ICount  =  1  To  INumUAVs 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  uavNames . Item (ICount) 

Next  ICount 

Set  missionNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet (rMissions,  missionNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (missions) .  User 

chose  to  cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumMissions  =  missionNames . Count 

For  ICount  =  1  To  INumMissions 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  missionNames . I tern ( ICount ) 

Next  ICount 

Set  payloadNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet (rPayloads,  payloadNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (payloads) .  User 

chose  to  cancel:  no  further  output." 
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GoTo  subExit 
End  If 

INumPayloads  =  pa yloadNames . Count 

For  ICount  =  1  To  INumPayloads 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  payloadNames . Item ( ICount) 

Next  ICount 

Set  classNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet ( rClasses ,  classNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (classes) .  User 
chose  to  cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumClasses  =  classNames . Count 

For  ICount  =  1  To  INumPayloads 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  payloadNames . Item ( ICount) 

Next  ICount 

Set  gcsNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet (rGCSs,  gcsNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (GCS) .  User  chose  to 
cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumGCSs  =  gcsNames . Count 

For  ICount  =  1  To  INumGCSs 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  gcsNames . Item (ICount) 

Next  ICount 

Set  IrsNames  =  New  StringDictionaryClass 

If  Def ineGroundedSet (rLRSs,  IrsNames)  =  False  Then 

Print  #logFN,  "Duplicate  element  in  grounded  set  (LRS) .  User  chose  to 
cancel:  no  further  output." 

GoTo  subExit 

End  If 

INumLRSs  =  IrsNames . Count 
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For  ICount  =  1  To  INumLRSs 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  IrsNames . Item (ICount) 
Next  ICount 

Set  periodNames  =  New  StringDictionaryClass 
per iodNames . Size  =  3  *  lPlanPeriods 
For  ICount  =  1  To  lPlanPeriods 

periodNames .Add  "t"  &  Format (ICount) 

Next  ICount 

For  ICount  =  lFirstPeriod  To  lLastPeriod 

Print  #logFN,  Format ( ICount )  &  "  :  "  &  periodNames . Item (ICount) 
Next  ICount 


t 


t 


{GIVEN}  PARAMETERS 

num_uavs ( 1 , v,  t )  number  of  uavs  available  at  lrs  1  in  tier  v  in  period 


num_mods (p, t)  number  of  payload  modules  available  of  ID  p  in  period 


gcs_cap (g, t, h) 
lrs_cap (1 , t) 
length (m) 
val (m) 
k (1,  s, t) 


max  num  UAVs  on  gcs  g  in  period  t  on  mission  type  h 
max  num  UAVs  at  LRS  1  in  period  t 
maximum  periods  of  coverage  of  mission  m 
value  of  mission  m 

indicator  if  lrs  1  used  by  schedule  s  in  period  t 


Print  #uasDataFN,  "  PARAMETERS" 

'  num_uavs  ( 1 ,  v,  t )  number  of  uavs  available  at  lrs  1  in  tier  v  in 
period  t 

Print  #uasDataFN,  "  num_uavs ( 1 , v,  t)  /" 

ReDim  INumUAVinPer iod ( 1  To  IrsNames . Count,  1  To  tierNames . Count, 
lFirstPeriod  To  lLastPeriod)  As  Long 

For  IPeriod  =  lFirstPeriod  To  lLastPeriod 

For  1LRS  =  1  To  IrsNames . Count 

For  ITier  =  1  To  tierNames . Count 

INumUAVinPeriod (1LRS,  ITier,  IPeriod)  =  0 

Next  ITier 

Next  1LRS 

For  iRow  =  1  To  rUAVs . Rows . Count 

If  rUAVs (iRow,  COL_UAV_FIRST_AVAIL)  <=  IPeriod  And  rUAVs (iRow, 
COL_UAV_LAST_AVAIL)  >=  IPeriod  Then 
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1LRS  =  IrsNames .Lookup (rUAVs (iRow,  COL_UAV_LRS) ) 

ITier  =  tierNames . Lookup ( rUAVs ( iRow,  COL_UAV_TIER) ) 

If  1LRS  =  0  Then 

MsgBox  "Error:  invalid  LRS  name  in  UAV  list" 
wsUAV. Activate 

rUAVs (iRow,  COLJJAVJLRS ). Select 
GoTo  subExit 
Elself  ITier  =  0  Then 

MsgBox  "Error:  invalid  tier  name  in  UAV  list" 
wsUAV. Activate 

rUAVs (iRow,  COL _UAV_TIER) . Select 
GoTo  subExit 
Else 

INumUAVinPeriod ( 1LRS ,  ITier,  IPeriod)  =  INumUAVinPeriod (1LRS, 
ITier,  IPeriod)  +  1 

End  If 

End  If 

Next  iRow 

Next  IPeriod 

For  1LRS  =  1  To  IrsNames . Count 

For  ITier  =  1  To  INumTiers 

For  IPeriod  =  lFirstPeriod  To  lLastPeriod 

If  INumUAVinPeriod ( 1LRS ,  ITier,  IPeriod)  >=  1  Then 

Print  #uasDataFN,  "  1_"  &  IrsNames . Item (1LRS)  &  ",v_"  & 

tierNames . Item ( ITier )  &  ",t_"  &  IPeriod  &  "  "  & 

INumUAVinPeriod (1LRS,  ITier,  IPeriod) 

End  If 

Next  IPeriod 
Next  ITier 
Next  1LRS 

Print  #uasDataFN,  "  /" 

ReDim  numMods (classNames . Count,  tierNames . Count )  As  Long 

For  IPayload  =  1  To  r Pa yloads . Rows . Count 

iRow  =  classNames . Lookup ( rPayloads ( IPayload,  COL_PAYLOAD_CLASS) ) 

iCol  =  tierNames . Lookup (rPayloads (IPayload,  COL_PAYLOAD_TIER) ) 

numMods (iRow,  iCol)  =  numMods ( iRow,  iCol)  +  rPayloads ( IPayload, 
COL_PAYLOAD_QUANT I T Y ) 

Next  IPayload 
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'  num_mods ( c, v)  number  of  payload  modules  available  of  class  c  for 

tier  v 

Print  #uasDataFN,  "  num_mods ( c, v)  /" 

For  lClass  =  1  To  classNames . Count 
For  ITier  =  1  To  tierNames . Count 

Print  #uasDataFN,  "c_"  &  classNames . I tern ( lClass )  &  ",v_"  & 

tierNames . Item ( ITier )  &  "  "  &  numMods (lClass,  ITier) 

Next  ITier 

Next  lClass 

Print  #uasDataFN,  "  /" 

'  gcs_cap (g, t, h)  max  num  UAVs  on  gcs  g  in  period  t  on  mission  type  h 

Print  #uasDataFN,  "  gcs_cap (g, t, h)  /" 

For  1GCS  =  1  To  INumGCSs 

For  IPeriod  =  lFirstPeriod  To  lLastPeriod 

Print  #uasDataFN,  "  g_"  &  gcsNames . Item (1GCS)  &  " . t_"  &  IPeriod 

&  " . h_active  "  &  rGCSs(lGCS,  COL_GCS_ACTIVE) 

Print  #uasDataFN,  "  g_"  &  gcsNames . Item ( 1GCS )  &  " . t_"  &  IPeriod 

&  " . h_passive  "  &  rGCSs  (1GCS,  C0LJ3CS  JPASSIVE) 

Next  IPeriod 

Next  1GCS 

Print  #uasDataFN,  "  /" 

'  lrs_cap(l,t)  max  num  UAVs  at  LRS  1  in  period  t 

Print  #uasDataFN,  "  lrs_cap(l,t)  /" 

For  1LRS  =  1  To  INumLRSs 

For  IPeriod  =  lFirstPeriod  To  lLastPeriod 

Print  #uasDataFN,  "  1_"  &  IrsNames . I tern ( 1LRS )  &  " . t_"  &  IPeriod 

&  "  "  &  rLRSs ( 1LRS ,  COL _LRS_ACTIVE) 

Next  IPeriod 

Next  1LRS 

Print  #uasDataFN,  "  /" 

'  length (m)  maximum  periods  of  coverage  of  mission  m 

Print  #uasDataFN,  "  length (m)  /" 

For  IMission  =  1  To  INumMissions 

Print  #uasDataFN,  "  m_"  &  missionNames .  Item  ( IMission)  &  "  "  & 

rMissions (IMission,  COL_MISSION_PERIODS ) 

Next  IMission 

Print  #uasDataFN,  "  /" 
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val (m) 


value  of  mission  m 


Print  #uasDataFN,  "  val (m)  /" 

For  IMission  =  1  To  INumMissions 

Print  #uasDataFN,  "  m_"  &  missionNames .  Item  ( IMission)  &  "  "  & 

rMissions (IMission,  COL_MISSION_VALUE) ~ 

Next  IMission 

Print  #uasDataFN,  "  /" 


'  Calculate  distances  between  GCS-Mission  pairs,  LRS-mission  pairs,  and 
mission-mission  pairs 

ReDim  dDistGM (INumGCSs,  INumMissions)  As  Double 
ReDim  dDistLM (INumLRSs,  INumMissions)  As  Double 
ReDim  dDistMM (INumMissions,  INumMissions)  As  Double 

Dim  i  As  Long,  j  As  Long 


For  1GCS  =  1  To  rGCSs . Rows . Count 

i  =  gcsNames . Lookup (rGCSs (1GCS,  COL_GCS_NAME) ) 

For  IMission  =  1  To  rMissions .Rows . Count 

j  =  missionNames . Lookup (rMissions (IMission,  COL_MISSION_NAME) ) 

dDistGM (i,  j)  =  SphericalDistance (rGCSs (1GCS,  COL_GCS_LAT) , 
rGCSs (1GCS,  COL_GCS  LON), 


rMissions (IMission, 

COL_MISSION_LAT)  ,  rMi ssions ( IMission,  COL_MISSION_LON) ) 

Next  IMission 


Next  1GCS 


For  1LRS  =  1  To  rLRSs . Rows . Count 

i  =  IrsNames. Lookup (rLRSs (1LRS,  COL_LRS_NAME) ) 

For  IMission  =  1  To  rMissions . Rows . Count 

j  =  missionNames . Lookup (rMissions (IMission,  COL_MISSION_NAME) ) 

dDistLM (i,  j)  =  SphericalDistance ( rLRSs ( 1LRS ,  COL_LRS_LAT) , 
rLRSs (1LRS,  COL_LRS  LON), 


rMissions (IMission, 

COL_MISSION_LAT)  ,  rMi ssions ( IMission,  COL_MISSION_LON) ) 

Next  IMission 


Next  1LRS 


For  IMission  =  1  To  rMissions .Rows . Count 

i  =  missionNames . Lookup (rMissions (IMission,  COL_MISSION_NAME) ) 
For  !Mission2  =  1  To  rMissions . Rows . Count 
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j  =  missionNames . Lookup (rMissions (lMission2 ,  COL_MISSION_NAME) ) 

dDistMM(i,  j)  =  SphericalDistance (rMissions (IMission, 

COL_MISSION_LAT) ,  rMissions (IMission,  COL_MISSION_LON) , 

rMissions (lMission2, 

COL_MISSION_LAT)  ,  rMissions (lMission2 ,  COL_MISSION_LON) ) 

Next  lMission2 
Next  IMission 

'Find  GCSs  within  range  of  each  mission.  If  none,  require  one  of  the 
AV_CLR  missions 

Dim  pFN  As  Integer,  hpFN  As  Integer 

Dim  numPreq  As  Long 

numPreq  =  0 

pFN  =  FreeFileO 

Open  "p.gms"  For  Output  As  pFN 

hpFN  =  FreeFileO 

Open  "hp.gms"  For  Output  As  hpFN 
Print  #pFN,  "  SET  P(m,mp)  /" 

Print  # hp FN ,  "  SET  HP (m)  /" 

For  IMission  =  1  To  missionNames . Count 
For  1GCS  =  1  To  gcsNames . Count 

If  dDistGM ( 1GCS ,  IMission)  <=  rGCSs (1GCS,  COL_GCS_LOS)  Then  GoTo 
nextMissionGCS 

Next  1GCS 

numPreq  =  numPreq  t  1 

Print  #hpFN,  "  m_"  &  missionNames . Item ( IMission) 

For  lMission2  =  1  To  rMissions . Rows . Count 

If  rMissions (lMission2,  COL_MISSION_CLASS ) 

wsGCS .Range ("LINK") .Value  Then 

Print  #pFN,  "  m_"  &  missionNames . Item (IMission)  &  ",m_"  & 

rMissions (lMission2,  COL_MISSION_NAME) 

End  If 

Next  lMission2 
nextMissionGCS : 

Next  IMission 
If  numPreq  =  0  Then 

Print  #pFN,  "  m_"  &  missionNames .  Item  (1 )  &  ".m_"  & 

missionNames . Item (1 ) 

Print  #hpFN,  "  m_"  &  missionNames . Item ( 1 ) 

End  If 

Print  #pFN,  "  /" 

Print  #pFN,  "  ; " 
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Close  pFN 

Print  # hp FN ,  "  /" 

Print  #  hp  FN ,  " 

Close  hpFN 

Dim  mgFN  As  Integer 
mgFN  =  FreeFileO 

Open  "mg.gms"  For  Output  As  mgFN 
Print  #mgFN,  "  SET  mg(m,g)  /" 

For  IMission  =  1  To  missionNames . Count 
For  iRow  =  1  To  rGCSs . Rows . Count 

1GCS  =  gcsNames . Lookup ( rGCSs ( iRow,  COL_GCS_NAME) ) 

If  dDistGM ( 1GCS ,  IMission)  <=  rGCSs (iRow,  COL_GCS_LOS)  Then 

Print  #mgFN,  "  m_"  &  missionNames . Item (IMission)  &  " . g_"  & 

gcsNames . Item (1GCS) 

End  If 

Next  iRow 

Next  IMission 

Print  #mgFN,  "  /;" 

Close  mgFN 

Dim  mhFN  As  Integer 
mhFN  =  FreeFileO 

Open  "mh.gms"  For  Output  As  mgFN 

ReDim  bClassActive ( classNames . Count )  As  Boolean 
For  iRow  =  1  To  rClasses . Rows . Count 

lClass  =  classNames . Lookup ( rClasses (iRow,  COL_MCLASS_NAME) ) 

If  CLng  (rClasses (iRow,  COL_MCLASS_ACTIVE) )  =  1  Then 
bClassActive (lClass)  =  True 
Else 

bClassActive (lClass)  =  False 
End  If 
Next  iRow 

Print  #mhFN,  "  SET  mh(m,h)  /" 

For  iRow  =  1  To  rMissions .Rows . Count 

lClass  =  classNames . Lookup ( rMissions ( iRow,  COL_MISSION_CLASS) ) 

If  bClassActive (lClass)  Then 
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Print  #mhFN,  "  m_"  &  rMissions (iRow,  COL_MISSION_NAME)  & 

" . h_active" 

Else 

Print  #mhFN,  "  m_"  &  rMissions (iRow,  COL_MISSION_NAME)  & 

" . h_passive" 

End  If 
Next  iRow 
Print  #mhFN,  "  /;" 

Close  mhFN 

Dim  gvFN  As  Integer 
gvFN  =  FreeFileO 

Open  "gv.gms"  For  Output  As  mgFN 
Print  #gvFN,  "  SET  GV(g,v)  /" 

For  iRow  =  1  To  rGCSs . Rows . Count 

1GCS  =  gcsNames . Lookup ( rGCSs ( iRow,  COL_GCS_NAME) ) 

ITier  =  tierNames . Lookup ( rGCSs ( iRow,  COL_GCS_TIER) ) 

Print  #gvFN,  "  g_"  &  gcsNames .  Item  ( 1GCS )  &  ",v_"  & 

tierNames . Item (ITier) 

Next  iRow 

Print  #gvFN,  "  /;" 

Close  gvFN 

Print  #uasDataFN,  "  r (m)  /" 

For  iRow  =  1  To  rMissions . Rows . Count 

IMission  =  missionNames . Lookup ( rMissions ( iRow,  COL_MISSION_NAME) ) 

If  Trim (CStr (rMissions (iRow,  COL_MISSION_REQUIRED) ) )  <>  ""  Then 

If  CLng  (rMissions (iRow,  COL_MISSION_REQUIRED)  )  > 

CLng (rMissions (iRow,  COL_MISSION_PERIODS) )  Then 

If  MsgBox ( "Required  coverage  exceeds  total  mission  time", 
vbOKCancel)  =  vbCancel  Then 

wsMissions .Activate 

rMissions (iRow,  COL_MISSION_REQUIRED) .Select 
Exit  Sub 
Else 

Print  #uasDataFN,  "  m_"  &  missionNames . Item (IMission)  &  "  " 

&  CLng (rMissions (iRow,  COL_MISSION_PERIODS) ) 

End  If 

Else 

Print  fuasDataFN,  "  m_"  &  missionNames . Item ( IMission)  &  "  "  & 

CLng (rMissions (iRow,  COL_MISSION_REQUIRED) ) 


End  If 
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Else 


Print  #uasDataFN,  "  m_"  &  missionNames . Item (IMission)  &  "  0" 

End  If 
Next  iRow 

Print  #uasDataFN,  "  /" 


Count  number  of  UAVs  in  each  tier  at  each  LRS 
ReDim  tierAtLRS (tierNames . Count,  IrsNames . Count)  As  Long 
For  1LRS  =  1  To  IrsNames . Count 
For  ITier  =  1  To  tierNames . Count 
tierAtLRS (ITier,  1LRS)  =  0 
Next  ITier 
Next  1LRS 


For  1UAV  =  1  To  rUAVs . Rows . Count 

i  =  tierNames .Lookup (rUAVs (1UAV,  COL_UAV_TIER) ) 
j  =  IrsNames. Lookup (rUAVs (1UAV,  COL_UAV_LRS) ) 
tierAtLRS (i,  j)  =  tierAtLRS (i,  j)  t  1 
Next  1UAV 


'  a  (m,  s,  t) 
period  t 

b  (v,  s, t) 

period  t 

f  (g,  s,  t,  h) 
mission  type  h 

'  k (1, s,  t) 

'  q(p,s,t) 


indicator  if  mission  m  covered  by  schedule  s  in 

indicator  if  UAV  in  Tier  v  required  by  schedule  s  in 

indicator  if  gcs  g  used  by  schedule  s  in  period  t  for 

indicator  if  lrs  1  used  by  schedule  s  in  period  t 
indicator  if  payload  p  used  by  schedule  s  in  period  t 


Dim  IMaxRoutes  As  Long 

IMaxRoutes  =  CLng (wsDashboard . Range ( "MAX_ROUTES " ) ) 

Dim  IMaxMissions  As  Long 

IMaxMissions  =  CLng  (wsDashboard . Range ( "MAX_MISSIONS ") )  t  1 

ReDim  a (IMaxRoutes,  wsDashboard. Range ("LAST_PERIOD") )  As  Long 
ReDim  b (IMaxRoutes)  As  Long 

ReDim  f_act (IMaxRoutes,  wsDashboard . Range ( "LAST_PERIOD" ) )  As  Long 
ReDim  f_pas (IMaxRoutes,  wsDashboard. Range ("LAST_PERIOD") )  As  Long 
ReDim  k (IMaxRoutes,  wsDashboard. Range ("LAST_PERIOD") )  As  Long 
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ReDim  q (IMaxRoutes,  classNames . Count)  As  Long 

ReDim  stack (IMaxMissions  +  1)  As  Long 
ReDim  onStack (INumMissions)  As  Long 
ReDim  dwell (IMaxMissions  +  1)  As  Long 

For  IMission  =  1  To  INumMissions 
onStack (IMission)  =  0 
Next  IMission 

Dim  top  As  Long 
top  =  0 

INumSchedules  =  0 
INumLRSSchedules  =  0 
bNextRoute  =  True 
1LRS  =  1 

Dim  s  As  Long,  t  As  Long 

Dim  INextMission  As  Long,  lCurrMission  As  Long 

ReDim  INextStack (IMaxMissions  +  1)  As  Long 

Dim  dRemTime  As  Double 

Dim  IRemPeriods  As  Long 

Dim  dTransitTime  As  Double 

Dim  ITransitPer iods  As  Long 

Dim  lcurrNumClasses  As  Long 

Dim  dTierSpeed  As  Double 

Dim  IThisLaunchPeriod  As  Long 

Dim  IThisRecoverPer iod  As  Long 

Dim  IMaxClasses  As  Long 

ReDim  ITransitPeriod (IMaxMissions  +  1)  As  Long 

ReDim  lldlePeriod (IMaxMissions  +  1)  As  Long 

ReDim  IMissionPeriod (IMaxMissions  +  1)  As  Long 

ReDim  IStartPer iod ( INumMissions )  As  Long 

ReDim  lFinishPeriod (INumMissions)  As  Long 

For  IMission  =  1  To  rMissions .Rows . Count 

IStart Period (missionNames . Lookup ( rMissions ( IMission, 
COL_MISSION_NAME) ) )  =  CLng ( rMis sions ( IMis sion ,  COL_MISSION_START_PERIOD) ) 

lFinishPeriod (missionNames . Lookup (rMissions (IMission, 
COL_MISSION_NAME) ) )  =  CLng ( rMis sions ( IMis sion ,  COL_MISSION_FINISH_PERIOD) ) 

Next  IMission 
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ReDim  dMissionValue (missionNames . Count)  As  Double 

For  IMission  =  1  To  rMissions . Rows . Count 

dMissionValue (missionNames . Lookup (rMissions (IMission, 

COL_MISSION_NAME) ) )  =  CDbl (rMissions (IMission,  COL_MISSION_VALUE) ) 

Next  IMission 

ReDim  lHasClass (INumClasses)  As  Long 
For  lClass  =  1  To  INumClasses 
lHasClass (lClass)  =  0 
Next  lClass 
lcurrNumClasses  =  0 
ReDim  reqClass (INumMissions) 

For  IMission  =  1  To  rMissions . Rows . Count 

reqClass  (missionNames . Lookup ( rMissions ( IMission,  COL_MISSION_NAME) ) ) 

classNames . Lookup ( rMissions ( IMission,  COL_MISSION_CLASS) ) 

Next  IMission 

ReDim  lClassTierAvail (classNames . Count,  tierNames . Count )  As  Long 
For  lClass  =  1  To  classNames . Count 
For  ITier  =  1  To  tierNames . Count 
lClassTierAvail  (lClass,  ITier)  =  0 
Next  ITier 
Next  lClass 

For  IPayload  =  1  To  rPayloads . Rows . Count 

lClass  =  classNames . Lookup ( rPayloads ( IPayload,  COL_PAYLOAD_CLASS) ) 

ITier  =  tierNames . Lookup ( rPayloads ( IPayload,  COL_PAYLOAD_TIER) ) 

lClassTierAvail (lClass,  ITier)  =  CLng ( rPayloads ( IPayload, 

COL_PAYLOAD_QUANT I T Y ) ) 

Next  IPayload 

Dim  aFN  As  Integer,  bFN  As  Integer,  fFN  As  Integer,  kFN  As  Integer,  qFN 
As  Integer 

aFN  =  FreeFileO 

Open  ThisWorkbook . path  &  "\a.gms"  For  Output  As  aFN 
bFN  =  FreeFileO 

Open  ThisWorkbook . path  &  "\b.gms"  For  Output  As  bFN 
kFN  =  FreeFileO 

Open  ThisWorkbook . path  &  "\k.gms"  For  Output  As  kFN 
qFN  =  FreeFileO 

Open  ThisWorkbook . path  &  "\q.gms"  For  Output  As  qFN 

s  =  0 

Print  #aFN,  "  SET  A(s,m,t)  /" 
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Print  #bFN,  "  SET  B(s,v,t)  /" 

Print  #qFN,  "  SET  Q(s,c,t)  /" 

Dim  tierLRScombos  As  Long 
tierLRScombos  =  0 
For  ITier  =  1  To  tierNames . Count 
For  1LRS  =  1  To  IrsNames . Count 
If  tierAtLRS (ITier,  1LRS)  Then 

tierLRScombos  =  tierLRScombos  +  1 
End  If 
Next  1LRS 
Next  ITier 

Dim  currCombo  As  Long 
currCombo  =  1 

For  1LRS  =  1  To  rLRSs . Rows . Count 

Print  #logFN,  "LRS  "  &  rLRSs (1LRS,  COL_LRS_NAME )  &  & 

IrsNames. Lookup (rLRSs (1LRS,  COL_LRS_NAME) ) 

For  ITier  =  1  To  rTiers . Rows . Count 

Print  #logFN,  "  tier  "  &  rTiers  ( ITier ,  COL_TIER_NAME)  &  & 

tierNames . Lookup (rTiers (ITier,  COL_TIER_NAME) ) 

If  tierAtLRS (tierNames . Lookup ( rTiers  ( ITier ,  COL_TIER_NAME) ) , 

IrsNames. Lookup (rLRSs (1LRS,  COL_LRS_NAME) ) )  >  0  Then 

dTierSpeed  =  rTiers  (ITier,  COL_TIER_SPEED) 

dRemTime  =  rTiers (ITier,  COL_TIER_ENDURANCE)  '*  rTiers (ITier, 

C0L_TIERJ3PEED) 

IRemPeriods  =  Ceiling (dRemTime  *  IPer iodsPerHour ) 

IMaxClasses  =  rTiers  ( ITier,  COL_TIER_CAPACITY) 

Print  #logFN,  "  tierAtLRS?  yes  lMaxClasses="  &  IMaxClasses 
For  lClass  =  1  To  INumClasses 
lHasClass (lClass)  =  0 
Next  lClass 
lcurrNumClasses  =  0 

top  =  1 

stack  (top)  =  EVENT_LR  ' L/R  event  at  1LRS 
INextStack (top)  =  1 
Do 

lCurrMission  =  stack (top) 

INextMission  =  INextStack (top) 

'Find  next  mission  to  put  on  the  stack 
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Do  While  INextMission  <=  INumMissions 
If  onStack ( INextMission)  >  0  Then 
'Print  #logFN,  "  onstack" 

GoTo  skipMission 
End  If 

If  dDistLM (1LRS,  INextMission)  >  rTiers (ITier, 

COL_TIER_RADIUS )  Then 

'Print  #logFN,  "Out  of  range" 

GoTo  skipMission 
End  If 

If  lClassTierAvail (reqClass (INextMission) ,  ITier)  =  0  Then 
'Print  #logFN,  "No  Payloads  available  for  this  Tier" 

GoTo  skipMission 
End  If 

If  lCurrMission  =  EVENT_LR  Then 

dTransitTime  =  dDistLM ( 1LRS ,  INextMission)  /  dTierSpeed 
Else 

dTransitTime  =  dDistMM (lCurrMission,  INextMission)  / 

dTierSpeed 

End  If 

ITransitPeriods  =  Ceiling (dTransitTime  *  IPeriodsPerHour ) 

If  lCurrMission  <=  0  Then 

If  lFinishPeriod (INextMission)  <= 

IStartPeriod (lCurrMission)  +  ITransitPeriods  Then 

'Print  #logFN,  "  time  warp" 

GoTo  skipMission 

End  If 

Elself  dMissionValue ( INextMission )  <— 

dMissionValue (lCurrMission)  Then 

If  lFinishPeriod (INextMission)  <= 

lFinishPeriod (lCurrMission)  +  ITransitPeriods  Then 

'Print  #logFN,  "  time  warp" 

GoTo  skipMission 

End  If 

Else 

If  IStartPeriod ( INextMission)  <=  IStartPeriod (lCurrMission) 
+  ITransitPeriods  Then 

'Print  #logFN,  "  time-value  warp" 

GoTo  skipMission 
End  If 
End  If 
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If  ITransitPeriods  +  Ceiling  (IPeriodsPerHour  *  (dDistLM (1LRS, 
INextMission)  /  rTier s ( ITier ,  COL_TIER_SPEED) )  )  >  IRemPeriods  Then 

'Print  #logFN,  "  out  of  time" 

GoTo  skipMission 

End  If 

If  ( lcurrNumClasses  =  IMaxClasses  And 

lHasClass (reqClass ( INextMission) )  =  0)  Then 

'Print  #logFN,  "  no  payload  space  left" 

GoTo  skipMission 

End  If 

GoTo  foundMission 

skipMission : 

INextMission  =  INextMission  +  1 
Loop 

foundMission : 

If  INextMission  >  INumMissions  Or  top  >=  IMaxMissions  Then  'out 
of  missions:  pop  stack 

If  lCurrMission  >  0  Then 

lHasClass (reqClass (lCurrMission) )  = 

lHasClass (reqClass ( lCurrMission ) )  -  1 

If  lHasClass (reqClass (lCurrMission) )  =  0  Then 

lcurrNumClasses  =  lcurrNumClasses  -  1 

End  If 

onStack (lCurrMission)  =  0 
End  If 

top  =  top  -  1 
GoTo  nextMission 
End  If 

If  lCurrMission  =  EVENT_LR  Then 

ITransitPeriod ( top)  =  ITransitPeriods 

IThisLaunchPeriod  =  IStartPeriod ( INextMission) 

ITransitPeriods  -  1 

If  IThisLaunchPeriod  <  1  Then 
IThisLaunchPeriod  =  1 
End  If 

IRemPeriods  =  Ceiling (dRemTime  *  IPeriodsPerHour)  + 
IThisLaunchPeriod  -  1 

ITransitPeriod ( top)  =  IThisLaunchPeriod  +  ITransitPeriods 
lldlePeriod (top)  =  IThisLaunchPeriod  -  1 
Else 

IMissionPeriod (top)  =  lFinishPeriod (lCurrMission) 
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And 


If  dMissionValue ( lCurrMission)  <  dMissionValue (INextMission) 


IStartPeriod (INextMission)  -  ITransitPeriods  1  < 

lFinishPeriod (lCurrMission)  Then 

IMissionPeriod (top)  =  IStartPeriod ( INextMission ) 

ITransitPeriods  -  1 

End  If 

ITransitPeriod ( top)  =  IMissionPeriod (top)  +  ITransitPeriods 
lldlePeriod (top)  =  IStartPeriod (INextMission)  -  1 
End  If 

INextStack ( top)  =  INextMission  +  1 
top  =  top  +  1 
stack (top)  =  INextMission 
onStack (INextMission)  =  top 

IMissionPeriod (top)  =  lFinishPeriod (INextMission) 

ITransitPeriod ( top)  =  lFinishPeriod (top)  + 

Ceiling ( (dDistLM ( 1LRS ,  INextMission)  /  rTier s ( ITier ,  COL_TIER_SPEED) )  * 
1 Pe r i ods Per Hour ) 

INextStack (top)  =  1 

If  lHasClass (reqClass (INextMission) )  =  0  Then 
lcurrNumClasses  =  lcurrNumClasses  +  1 
End  If 

lHasClass (reqClass (INextMission) )  = 

lHasClass (reqClass ( INextMission) )  +  1 

'New  feasible  schedule:  now  copy  stack  info  into  appropriate 

arrays 

s  =  s  +  1 
t  =  1 

'calculate  start  period  of  mission  sequence:  find  start  time  of 

first  mission 

'  in  stack (2),  figure  out  latest  departure  period  to  arrive  at 

first  mission  (if 

'  earlier  than  period  1,  then  1/r  in  period  1,  transit,  and  do 
as  much  of  mission 

'  as  possible 

Print  #logFN,  "dTransitTime  "  &  dTransitTime  &  " 

ITransitPeriods  "  &  ITransitPeriods 

Print  #logFN,  "s:"  &  CStr(s)  &  " 

For  iRow  =  1  To  top 

Print  #logFN,  stack  (iRow)  & 

Next  iRow 

Print  #logFN,  "0" 
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Print  #logFN,  "m: "  &  CStr(s)  &  " 

For  iRow  =  1  To  top 

Print  #logFN,  IMissionPeriod (iRow)  & 

Next  iRow 
Print  #logFN,  "0" 

Print  #logFN,  "t:"  &  CStr(s)  &  " 

For  iRow  =  1  To  top 

Print  #logFN,  ITransitPeriod ( iRow)  & 

Next  iRow 
Print  #logFN,  "0" 

Print  #logFN,  "i:"  &  CStr(s)  &  " 

For  iRow  =  1  To  top 

Print  #logFN,  lldlePeriod (iRow)  & 

Next  iRow 
Print  #logFN,  "0" 

Do  While  t  <  IThisLaunchPeriod  And  t  <=  lLastPeriod 
a (s,  t)  =  EVENT_IDLE 
t  =  t  +  1 
Loop 

a  (s,  t)  =  EVENT_LR 
f_act (s,  t)  =1 
k ( s ,  t)  =  1LRS 
t  =  t  t  1 

Do  While  t  <=  ITransitPeriod ( 1 )  And  t  <=  lLastPeriod 
a (s,  t)  =  EVENT_TRANSIT 
f_pas(s,  t)  =  1 
t  =  t  +  1 
Loop 

For  iRow  =  2  To  top 

lCurrMission  =  stack (iRow) 

If  iRow  <  top  Then 

INextMission  =  stack (iRow  +  1) 

Else 

INextMission  =  EVENT_LR  'L/R  event 
End  If 

'  XXX  Need  upper  bounds:  time  horizon,  flight  time,  etc... 
Do  While  t  <=  IMissionPeriod (iRow)  And  t  <=  lLastPeriod 
' fill  in  arrays 
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a(s,  t)  =  lCurrMission 
f_act  (s,  t)  =1 
t  =  t  +  1 
Loop 

Do  While  t  <=  ITransitPer iod ( iRow)  And  t  <=  lLastPeriod 
a (s,  t)  =  EVENT_TRANSIT 
f_pas(s,  t)  =1 
t  =  t  +  1 
Loop 

Do  While  t  <=  lldlePeriod (iRow)  And  t  <=  lLastPeriod 
a (s,  t)  =  EVENT_IDLE 
f_pas(s,  t)  =1 
t  =  t  +  1 
Loop 

Next  iRow 

If  t  <=  lLastPeriod  Then 
a (s,  t)  =  EVENT_LR 
IThisRecoverPeriod  =  t 
f_act  (s,  t)  =1 
k  ( s ,  t)  =  1LRS 
b ( s )  =  ITier 
End  If 

For  iRow  =  1  To  lLastPeriod 

Print  #logFN,  &  a(s,  iRow); 

Next  iRow 

Print  #logFN,  "" 

For  iRow  =  2  To  top 

For  IPeriod  =  1  To  lLastPeriod 

If  a(s,  IPeriod)  =  stack (iRow)  Then 

Print  #aFN,  "  s_"  &  CStr(s)  &  ",m_"  & 

missionNames . Item (stack (iRow) )  &  " . t_"  &  CStr ( IPeriod) 

End  If 

Next  IPeriod 
Next  iRow 

For  IPeriod  =  1  To  lLastPeriod 

If  IPeriod  >=  IThisLaunchPeriod  And  IPeriod  <= 
IThisRecoverPeriod  Then 

Print  #bFN,  "  s_"  &  CStr(s)  &  ",v_"  & 

Trim (rTiers (ITier,  COL_TIER_NAME) )  &  " . t  "  &  CStr ( IPeriod) 
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End  If 


Next  IPeriod 

Print  #kFN,  "  k ( ' s_"  &  CStr(s)  &  "','1  "  &  Trim (rLRSs (1LRS, 

COL_LRS_NAME) )  &  " ' , ' t_"  &  CStr ( IThisLaunchPeriod)  &  '")=1;" 

Print  #kFN,  "  k ( ' s_"  &  CStr(s)  &  " ' , ' 1_"  &  Trim (rLRSs (1LRS, 

COL_LRS_NAME) )  &  " ' , ' t_"  &  CStr ( IThisRecoverPeriod)  &  '")=1;" 

For  lClass  =  1  To  classNames . Count 

If  lHasClass ( lClass )  >  0  Then 

q(s,  lClass)  =  1 

For  IPeriod  =  1  To  lLastPeriod 

If  IPeriod  >=  IThisLaunchPeriod  And  IPeriod  <= 
IThisRecoverPeriod  Then 

Print  #qFN,  "  s_"  &  CStr(s)  &  " . c_"  & 

classNames . Item (lClass)  &  ".t_"  &  CStr ( IPeriod) 

End  If 

Next  IPeriod 
End  If 
Next  lClass 

If  s  >=  IMaxRoutes  Then  GoTo  doneMissions 

If  s  Mod  1000  =  0  Then  wsDashboard. Range ( "ROUTESCELL" )  =  s 
If  s  >=  (currCombo  -  1)  *  IMaxRoutes  /  (tierLRScombos)  Then 

GoTo  nextTier 

nextMission : 

Loop  While  top  >  0 
End  If 
nextTier : 

currCombo  =  currCombo  t  1 
Next  ITier 
Next  1LRS 

doneMissions : 

wsDashboard . Range ( "ROUTESCELL" )  =  s 

Print  #uasSetsFN,  "  s  /s_l*s_"  &  CStr(s)  &  "  /" 

Print  #uasSetsFN,  "  ; " 

Print  #uasDataFN,  "  ; " 

Print  #aFN,  "  /" 

Print  #aFN,  "  ; " 

Print  #bFN,  "  /" 

Print  #bFN,  "  ; " 

Print  #qFN,  "  /" 

Print  #qFN,  "  ; " 
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subExit : 


Print  #logFN,  "Ending  output" 
Close  uasSetsFN 
Close  uasDataFN 
Close  aFN 
Close  bFN 
Close  kFN 
Close  qFN 
Close  logFN 
End  Sub 
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APPENDIX  C  GAMS  CODE 


$TITLE  UAS  Planner  08.06.02 


$ontext 

MCCDC/OAD  UAS  Mission  Planner 
Assigns  a  fleet  of  UAVs  to  a  list  of  missions  with  various 
requirements  and  time  windows  over  a  fixed  time  horizon  of  many 
time  periods. 

Solved  using  explicit  column  generation,  one  column  per  UAS  schedule. 

Sofftext 

$inlinecom  {  } 

FILE  UASP  log/UASP.log/ ; 

PUT  UASPlog ; 

PUT  TJASP  -  MCCDC/OAD  UAS  Planner  08.06.02’  //  ; 

PUTCLOSE  UASP  log  ; 

FILE  UASP_sta/UASP .  sta/  ; 

PUT  UASP_sta ; 

PUT  ’ERROR  -  Did  not  solve'  //  ; 

PUTCLOSE  UASP  sta ; 

FILE  VALUE_CSV/VALUE.csv/ ; 

PUT  VALUECSV  ; 

PUT  ’  ’  /  ;  {  scratch  legacy  file  contents  } 

PUTCLOSE  VALUE  CSV  ; 

FILE  MIS  SION  SC  S  V/MIS  SION  S .  csv/ ; 

PUT  MISSIONSCSV  ; 

PUT  '  ’  /  ;  {  scratch  legacy  file  contents  } 

PUTCLOSE  MISSIONS  CSV ; 
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FILE  SCHEDULES_CSV/SCHEDULES.csv/ ; 
PUT  SCHEDULESCSV  ; 

PUT  ’ '  /  ;  {  scratch  legacy  file  contents  } 

PUTCLOSE  SCHEDULES  CSV 

FILE  DUALS_CSV/DUALS.csv/ ; 

PUT  DUALS  CSV ; 

PUT  ’ '  /  ;  {  scratch  legacy  file  contents  } 

PUTCLOSE  DUALS  CSV 


PUT  UASPJog ; 

UASP_log.ap=l  ;  {  re-open  log  file,  appending  to  prior  contents  } 

OPTIONS 
op  ter  =0.1 
liinrow  =  0 

limcol  =  0 
SEED  =  1234 

ITERLIM  =  10000000 
LP  =  CPLEX 
RMIP  =  CPLEX 

MIP  =  CPLEX 

SONTEXT 

{STATIC,  GROUNDED}  SETS 
v  Tier  levels 

h  mission  types  {active,  passive} 

g  Ground  control  stations 

1  Launch  and  recovery  sites 

m  Missions  {missions  that  ships  can  be  assigned} 

c  Payload  classes 
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t 


Time  periods 

UAV  emplyment  schedules 


s 

$OFFTEXT 

$INCLUDE  uasSets.gms 

alias(m,mp); 

$ONTEXT 

{GIVEN}  PARAMETERS 

num_uavs(l,v,t)  number  of  uavs  available  in  tier  v  in  period  t 
num_mods(c,t)  number  of  payload  modules  available  of  class  c  in  period  t 
gcs_cap(g,t,h)  max  num  UAVs  on  gcs  g  in  period  t  on  mission  type  h 

lrs_cap(l,t)  max  num  UAVs  at  LRS  1  in  period  t 

length(m)  maximum  periods  of  coverage  of  mission  m 

val(m)  value  of  mission  m 

r(m)  required  periods  of  coverage  for  mission  m 

k(s,l,t)  indicator  if  lrs  1  used  by  schedule  s  in  period  t 

SOFFTEXT 
SOFFLISTING 
$INCLUDE  uasData.gms 
PARAMETER 
k(s,l,t) 

$  INCLUDE  k.gms 

SONTEXT 
{DYNAMIC}  SETS 

A(s,m,t)  schedule  s  covers  mission  m  in  period  t 
B(s,v,t)  schedule  s  uses  tier  v  UAV  in  period  t 
Q(s,c,t)  schedule  s  requires  payload  class  c  in  period  t 
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P(m,mp)  mission  prerequisite  fixed  mission  m  requires  at  least  one 
mission  mp  in  same  period 

E(m,mp)  exclusive  missions  m  and  mp  cannot  be  accomplished  in  same 
period 

{Notional  data:  will  be  replaced  with  INCLUDE  statements  and  accompanying 
data  files} 

$OFFTEXT 

alias(m,mp); 

alias(v,vp); 

alias(t,tp); 

$INCLUDE  a.gms 

$  INCLUDE  b.gms 

$  INCLUDE  q.gms 

$  INCLUDE  hp.gms 

$  INCLUDE  p.gms 

$INCLUDE  mg.gms 

$  INCLUDE  mh.gms 

$  INCLUDE  gv.gms 

$ONLISTING 

SET  sl(s,l); 

LOOP((s,l)$(SUM(t,k(s,l,t))>eps), 
sl(s,l)  =  yes; 

); 


SCALARS  errors, warnings  ; 

errors=0  ; 
warnings=0; 

{Loops  for  data  checking  go  here} 

IF(errors>0, 

PUT  /  ’inconsistencies  found  in  active  and  value  data:  terrors:  10:0  / ; 
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PUT  '(should  I  have  shut  you  down  here?'  /  ; 

); 


VARIABLE 

OBJECTIVE 

BINARY  VARIABLES 
X(s)  Schedule  s  flown  by  some  UAV 

Y(m,v,t)  Mission  m  covered  in  period  t 
W(m,v,g,t)  Mission  m  handled  by  GCS  g  in  period  t 

9 

INTEGER  VARIABLES 

D(m)  Total  number  of  dwell  periods  on  mission  m 
MODS(c,v,l) 

9 

EQUATIONS 

eqnObj  {RO}  Objective  function  measures  total  value  of  mission 

accomplishment 

UAVLimit(l,v,t)  {R1 }  limit  on  UAVs  flying  by  tier  and  period 

ModuleLimit(c,v,l,t)  {R2}  limit  on  payload  modules  flying  by  payload  and 

period 

TotalMods(c,v)  {R2.5} 

GCSLimit(g,t,h)  {R3 }  limit  on  GCS  use  by  gcs  period  and  type 

LRSLimit(l,t)  {R4}  limit  on  LRS  use  by  LRS  and  period 

CoverageLimit(m,v,t)  {R5}  mission  and  period  coverage  control 
DwellLimit(m)  {R6}  limit  on  total  dwell  time  by  mission 

Prerequisite(m,v,t)  {R7}  prerequisite  controls  by  period 
Mission_Control(m,v,t) 

*  Exclusive(m,mp,t)  {R8}  mutually  exclusive  missions  by  period 
RequiredDwell(m)  {R9}  required  missions  must  be  completely  covered 
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eqnObj..  {RO} 

OBJECTIVE  =e=  sum(m,  val(m)*D(m)) 


UAVLimit(l,v,t)..  {R1 } 

SUM(s$(B(s,v,t)  and  sl(s,l)),  X(s))  =1=  num_uavs(l,v,t) 

ModuleLimit(c,v,l,t)..  {R2} 

SUM(s$(sl(s,l)  and  B(s,v,t)  and  Q(s,c,t)),  X(s))  =1=  MODS(c,v,l) 

9 

TotalMods(c,v).. 

SUM(1,  MODS(c,v,l))  =1=  num_mods(c,v) 

9 

GCSLimit(g,t,h)..  {R3} 

SUM((m,v)$(MG(m,g)  and  MH(m,h)  and  GV(g,v)),  W(m,v,g,t)) 

gcs_cap(g,t,h) 

9 

LRSLimit(l,t)..  {R4} 

SUM(s$(k(s,l,t)>eps),  X(s))  =1=  lrs_cap(l,t) 

9 

CoverageLimit(m,v,t)..  {R5 } 

Y(m,v,t)  =1=  SUM(s$(A(s,m,t)  and  B(s,v,t)),  X(s)) 

9 

DwellLimit(m)..  {R6} 

D(m)  =1=  SUM((v,t),  Y(m,v,t)) 

9 

Prerequisite(m,v,t)$HP(m)..  {R7} 

Y(m,v,t)  =1=  SUM((mp,vp)$P(m,mp),Y(mp,vp,t)) 

9 

Mission_Control(m,v,t)$(not  HP(m)).. 

Y(m,v,t)  =1=  SUM(g$(MG(m,g)  and  GV(g,v)),  W(m,v,g,t)) 

9 

*  Exclusive(m,mp,t)$E(m,mp)..  {R8} 
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Y(m,t)  +  Y(mp,t)  =1=  1 


* 

*  . 

? 

RequiredDwell(m)$(r(m)>eps). .  {R9 } 

D(m)  =g=  r(m) 

LOOP(m, 

D.up(m)  =  length(m);  {RIO} 

); 


MODEL  UAS_PLANNER  / 
eqnObj 
UAVLimit 
ModuleLimit 
TotalMods 
GCSLimit 
LRSLimit 
CoverageLimit 
DwellLimit 
Prerequisite 
*  Exclusive 
MissionControl 
RequiredDwell 

/ 


IF(errors=0, 

PUT  'solve  UAS_PLANNER’  /  ; 

*  UAS_PLANNER.reslim=60.  ; 

*  U AS_PL ANNER.  optcr=0 . 0  5  ; 

PUT  ’  reslim=',UAS_PLANNER.reslim:5:0,' 

optcr=',UAS_PLANNER.optcr:5:2  ; 

PUT  ’  optfde=',UAS_PLANNER.optfile:2:0  / ; 


seconds, 
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SOLVE  UAS  PLANNER  USING  MIP  MAXIMIZING  OBJECTIVE; 

PUT  'PLANNER  solver  and  model 

',UAS_PLANNER.solvestat:4:0,UAS_PLANNER.modelstat:4:0  / ; 

IF(  UAS_PLANNER.modelstat<>  1  and  UAS_PLANNER.modelstat<>8, 

errors=l  ; 

ELSE 

PUT  /  'Objective:  ',OBJECTIVE.1:10:4  /; 

PUT  /  'Module  Deployment'  /  ; 

LOOP((c,v,l)$(MODS.L(c,v,l)>eps), 

PUT  ’  ',c.tl:  10,'  ',v.tl:8,'  ',l.tl:  10,':  ',MODS.L(c,v,l):4:0  / ; 

); 

PUT  /  'UAV  flight  schedules..'  // ; 

LOOP(s$(X.l(s)>eps), 

PUT’  ',s.tl:  10,'  '; 

LOOP(v$SUM(t$B(s,v,t),  1), 

PUT  ’  ',v.tl:8  ; 

); 

PUT’  ',X.l(s):4:l  / ; 

LOOP(m$(SUM(t$  A(s,m,t),  1  )>0), 

PUT’  ',m.tl:10 ; 

LOOP(t$A(s,m,t), 

PUT  ’  ',t.tl:5  ; 

); 

PUT/; 

); 

); 

PUT  /  'GCS  Usage’  / ; 

LOOP(g$(SUM((m,v,t)$(MG(m,g)  and  GV(g,v)),W.l(m,v,g,t))>eps), 
PUT’  ',g.tl:8  / ; 

LOOP((m,v)$(MG(m,g)  and  GV(g,v)  and  (SUM(t,W.l(m,v,g,t))>eps)), 
PUT  ’  ',m.tl:8,'  ’,(SUM(t,W.l(m,v,g,t)))  ; 

LOOP(t$(W.l(m,v,g,t)>eps), 


status 
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PUT  ’  ',t.tl:8  ; 

); 

PUT/; 

); 

); 

PUT  /  ’MISSION  COVERAGE’  /  ; 

PUT  ’  Mission  Dwell  Time’  /  ; 

LOOP(m$(D.l(m)>0), 

PUT’  ’,m.tl:8,’  ’,D.l(m):5:l  /  ; 

LOOP(t$(SUM(v,Y.l(m,v,t))>eps), 

PUT  ’  ',t.tl:4,'  ',(SUM(v,Y.l(m,v,t))):4: 1  /  ; 

); 

); 

PUTCLOSE  UASP  LOG; 

PUT  MISSIONSCSV  ; 

PUT  ',„Total,Total, Value, Total’  /; 

PUT  ’, Start, End,Time  Periods, Time  Periods, Per  Time, Value’  /; 
PUT 

’Mission,Time, Time, Requested, Assigned, Period, Achieved,', OB JECTIVE.L  /; 
LOOP(m, 

PUT  m.tl,'„„',D.l(m),',',val(m),',',(D.l(m)*val(m))  /; 

); 

PUTCLOSE  MISSIONS  CSV; 

{Rest  of  report  writer  goes  here} 

); 


IF(errors>0, 

PUT  '  <***  errors  aborting  this  run:', errors:  10:0  /  ; 

); 

IF(warnings>0, 
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PUT/’  <+++  warnings  generated  by  this  run:  warnings:  10:0  / 
PUT  ’  please  find  out  what  these  indicate'  /  ; 

); 

PUTCLOSE  UASP  log  ; 


IF(errors=0, 


{Build  CSV  Output  Files} 


PUT  UASPSTA; 

PUT  'Optimal  solution  found.  Successful  completion.’  /; 
PUTCLOSE  UASP  STA; 

); 
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